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ABSTRACT 

The Planck mission is the most sensitive all-sky CMB experiment currently planned. The 
High Frequency Instrument {HFI) will be especially suited to observe clusters of galaxies by 
their thermal Sunyaev-Zel'dovich (SZ) effect. In order to assess Planck's, SZ-capabilities in 
the presence of spurious signals, a simulation is presented that combines maps of the thermal 
and kinetic SZ-effects with a realisation of the cosmic microwave background (CMB), in 
addition to Galactic foregrounds (synchrotron emission, free-free emission, thermal emission 
from dust, CO-line radiation) as well as the sub-millimetric emission from celestial bodies of 
our Solar system. Additionally, observational issues such as the finite angular resolution and 
spatially non-uniform instrumental noise of Planck's, sky maps are taken into account, yielding 
a set of all-sky flux maps, the auto-correlation and cross-correlation properties of which are 
examined in detail. In the second part of the paper, filtering schemes based on scale-adaptive 
and matched filtering are extended to spherical data sets, that enable the amplification of 
the weak SZ-signal in the presence of all contaminations stated above. The theory of scale- 
adaptive and matched filtering in the framework of spherical maps is developed, the resulting 
filter kernel shapes are discussed and their functionality is verified. 

Key words: galaxies: clusters: general, cosmology: cosmic microwave background, methods: 
numerical, space vehicles: Planck 



1 INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect jSunvaev &Zel'doviclJl97l 
ll98d:lRephaelill995irBirkinsha\\j|l999h is the most important ex- 
tragalactic source of secondary anisotropies in the CMB sky. The 
thermal SZ-effect is explained by the fact that CMB photons are put 
in thermal contact with electrons of the hot intra-cluster medium 
(ICM) by Compton-interactions which causes a transfer of energy 
from the ICM to the CMB. Because of the smallness of the Thomp- 
son cross-section and of the diluteness of the ICM this transfer of 
thermal energy is small. In the direction of a cluster, low-energetic 
photons with frequencies below v = 217 GHz are removed from 
the line-of-sight. At frequencies above v = 217 GHz CMB photons 
are scattered into the line-of-sight, causing a distinct modulation of 
the CMB surface brightness as a function of observing frequency, 
which enables the detection of clusters of galaxies in microwave 
data. 

In contrast, in the kinetic effect it is the peculiar motion of a 
cluster along the line of sight relative to the CMB frame that in- 
duces CMB surface brightness fluctuations. The peculiar motion of 
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the cluster causes the CMB to be anisotropic in the cluster frame. 
Due to this symmetry breaking of the scattering geometry, photons 
scattered into the line-of-sight are shifted in frequency, namely to 
higher frequencies, if the cluster is moving towards the observer. 

The P/anc/r-mission will be especially suited to detect SZ- 
clusters due to its sensitivity, its spectroscopic capabilities, sky cov- 
erage and spatial resolution. It is expected to yield a cluster cata- 
logue containing ^ 10** entries. Extensive literature exists on the 
topic, but so far the influence of foregrounds and details of Planck's 
instrumentation and data aquisition have not been thoroughly ad- 
dressed. In this work we aim at modelling the astrophysical and 
instrumental issues connected to the observation of SZ-clusters as 
exhaustively as possible: A simulation is presented that combines 
realistic maps of both SZ-efFects with a realisation of the CMB, 
with four different Galactic foreground components (thermal dust, 
free-free emission, synchrotron emission and emission from rota- 
tional transitions of CO molecules), with maps containing the sub- 
millimetric emission from planets and asteroids of the Solar sys- 
tem and with instrumental noise. Planck's frequency response and 
beam shapes are modelled conforming to the present knowledge 
of Planck's receivers and its optical system. In order to extract the 
SZ-cluster signal, filtering schemes based on matched and scale- 
adaptive filtering are extended to spherical data sets. 

In contrast to the recent work by iGeisbusch et alj i2004) . 
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our SZ-simulation does not rely on idealised scaling relations 
and futhermore, takes account of the cluster's morphological va- 
riety. The Galactic foregrounds are modelled in concordance with 
WMAP observations (see Bennett et al. 2003), which constitues an 
improvement over the simplifying assumptions made by Geisbiisch 
et al. In addition, instrumentation issues such as non-isotropic de- 
tector noise are properly incorporated into the simulation. The filter 
scheme employed in the paper by Geisbiisch et al. is the harmonic- 
space maximum entropy method introduced by Stolv arov et alj 
y003) which assumes approximate prior knowledge of the emis- 
sion component's power spectra. Its computational demand is much 
higher than matched and scale-adaptive filtering: In fact, the com- 
putations presented in this work can be run on a notebook-class 
computer. 

The paper is structured as follows: After a brief recapitulation 
of the SZ-effect in Sect.|5| the P/a«ffc-satellite and instrumental is- 
sues connected to observation of CMB anisotropics are decribed 
in Sect.0 The foreground emission components are introduced in 
Sect. 13] The steps in the simulation of flux maps for the various 
/'/ancA;-channels are described and their correlation properties are 
examined in Sect.|5| The theory of matched and scale-adaptive fil- 
tering is extended to spherical data sets and the resulting filter ker- 
nel shapes are in detail discussed in Sect.|6| A summary in Sect.|7| 
concludes the paper. 

Throughout the paper, the cosmological model assumed is the 
standard ACDM cosmology, which has recently been supported 
by observations of the WMAP satelhte i Spergel e tal. 2003). Pa- 
rameter values have been chosen as Om = 0.3, = 0.7, Hq = 
100/jkm S-' Mpc"' with h = 0.7, Ob = 0.04, = I and erg = 0.9. 



2 SUNYAEV-ZEL'DOVICH DEFINITIONS 

The Sunyaev-Zel'dovich effects are the most important extragalac- 
tic sources of secondary anisotropics in the CMB. Inverse Compton 
scattering of CMB photons with electrons of the ionised ICM gives 
rise to these effects and induce surface brightness fluctuations of 
the CMB sky, either because of the thermal motion of the ICM 
electrons (thermal SZ-effect) or because of the bulk motion of the 
cluster itself relative to the comoving CMB-frame along the line- 
of-sight (kinetic SZ-effect). 

The relative change AT/T in thermodynamic CMB temper- 
ature at position as a function of dimensionless frequency x = 
hv I {kgT cms) due to the thermal SZ-effect is given by: 
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where the amplitude y of the thermal SZ-effect is commonly known 
as the thermal Comptonisation parameter, that itself is defined as 
the line-of-sight integral of the temperature weighted thermal elec- 
tron density, m^, c, kg and crj denote electron mass, speed of light, 
Boltzmann's constant and the Thompson cross section, respec- 
tively. The kinetic SZ-effect arises due to the motion of the cluster 
parallel to the line of sight relative to the CMB-frame: 
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Here, v,. is the radial component of the cluster's velocity. The con- 
vention is such that v,. < 0, if the cluster is moving towards the ob- 
server. In this case, the CMB temperature is increased. In analogy. 



the quantity w is refered toas the kinetic Comptonisation. The SZ- 
observables are the line-of-sight Comptonisations integrated over 
the solid angle subtended by the cluster. The quantities if and "W 
are refered to as the integrated thermal and kinetic Comptonisa- 
tions, respectively: 
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Here, rfA(z) denotes the angular diameter distance of a cluster situ- 
ated at redshift z- 



3 SUBMILLIMETRIC OBSERVATIONS WITH Planck 

The P/ancfc-mission' - will perform a polarisation sensitive survey 
of the complete microwave sky in nine observing frequencies from 
the Lagrange point L2 in the Sun-Earth system. It will observe 
at angular resolutions of up to 5'.0 in the best channels and will 
achieve micro-Kelvin sensitivity relying on bolometric receivers 
(high frequency instrument HFI, described in Lamarre et al. 20(^ 
and on high elect ron mobility transistors (low frequency instrument 
LFI, see lVilla et al . 2003; BersanelhA Mandojesi 2000). The main 
characteristics are summarised in Table0 Planck's beam character- 
istics are given Sect. B.ll and the scanning strategy and the simula- 
tion of spatially non-uniform detector noise is outlined in Sect. l3.2l 



3.1 Beam shapes 

The beam shapes of Planck are well described by azimuthally sym- 
metric Gaussians b(0) = exp ( - | with cto = .^"^^ . The 

residuals from the ideal Gaussian shape (ellipticity, higher order 
distortions, diffraction rings, far-side lobes, pick-up of stray-light) 
are expected not to exceed the percent level and are neglected for 
the purpose of this work. TableQgives the angular resolution A6 in 
terms of FWHM of each Planck-channel for reference. 



3.2 Scanning strategy and noise-equivalent maps 

CMB observations by Planck will proceed in great circles fixed 
on the ecliptic poles. A single scan will start at the North ecliptic 
pole, will follow a meridian to the South ecliptic pole and back to 
the North ecliptic pole by following the antipodal meridian. Such 
a scan will last one minute and will be repeated sixty times. After 
that, the rotation axis will be shifted in a precessional motion for 
2^5 (approximately half a beam diameter) and the scan repeated. In 
this way, the entire sky is mapped once in 180 days. 

Fourier transform of the noise time series of Planck's receivers 
yields a noise power spectrum P(f) of the shape 
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i.e. the noise consists of two components: a power law component 
in frequency /, decribed by the spectral index a that assuming val- 
ues < Q- < 2 and a white noise component, smoothly joined at the 
frequency /knee- 

The /"-part of the noise spectrum originates from zero point 
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Planck channel 


1 


2 


3 


4 


5 


6 


1 


8 


9 


centre frequency vq 


30 GHz 


44 GHz 


70 GHz 


100 GHz 


143 GHz 


217 GHz 


353 GHz 


545 GHz 


857 GHz 


frequency window Ay 


3.0 GHz 


4.4 GHz 


7.0 GHz 


16.7 GHz 


23.8 GHz 


36.2 GHz 


58.8 GHz 


90.7 GHz 


142.8 GHz 


resolution A6» (FWHM) 


33^4 


26^8 


13fl 


9'.2 


TA 


5fO 


5fO 


5;o 


5fO 


noise level o"n 


1.01 mK 


0.49 mK 


0.29 mK 


5.67 mK 


4.89 mK 


6.05 mK 


6.80 mK 


3.08 mK 


4.49 mK 


thermal SZ-flux (Sy) 


-12.2 Jy 


-24.8 Jy 


-53.6 Jy 


-82.1 Jy 


-88.8 Jy 


-0.7 Jy 


146.0 Jy 


76.8 Jy 


5.4 Jy 


kinetic SZ-flux (5-^/) 


6.2 Jy 


13.1 Jy 


30.6 Jy 


55.0 Jy 


86.9 Jy 


llO.OJy 


69.1 Jy 


15.0 Jy 


0.5 Jy 


antenna temperature ATy 


-440 nK 


-417 nK 


-356 nK 


-267 nK 


-141 nK 


-0.5 nK 


38 «K 


8.4 nK 


0.2 nK 


antenna temperature AT'^y 


226 nK 


220 nK 


204 mK 


179 nK 


138 nK 


76 nK 


18 «K 


1.6 «K 


0.02 nK 



Table 1. Characteristics of Planck' & LF/-receivers (column 1-3) and ///"/-bolometers (column 4-9): centre frequency vq, frequency window Av as defined in 
eqns. 1271 and I28i . angular resolution A6 stated in FWHM, effective noise level ctn, fluxes (Sj/) and (Stv) generated by the respective Comptonisation of 
J/ = 'W = 1 arcmin" and the corresponding changes in antenna temperature ATy and AT^w- Due to Planck's, symmetric frequency response window, the 
thermal SZ-elfect does not vanish entirely at v = 217 GHz. 



drifts of the detector gain on large time scales. This power law com- 
ponent exhibits low-frequency variations that lead to the typical 
strip e pattern in s imulated Planck-ma.p& due to the scanning strat- 
egy (Maino et al. 1999). Algorithms for destriping the maps are 
a current research topic (for exa mple, the Mir age-algorithm pro- 
posed bv lYvon&Mavell j2004) . M APCimBA bvlDore et alJ <2()0ll) 
and the max-likelihood algorithm bv lNatoli et alJ but it can 

be expected that the destriping can be done very efficiently such that 
the remaining noise largely consists of uncorrelated pixel noise. 

In order to incorporate uncorrelated pixel noise into the sim- 
ulation, a set of maps has been construced, where at each pixel 
a number from a Gaussian distribution with width ctn has been 
drawn. For Planck'?, /fF/-receivers, the rms-fluctuations ctn in an- 
tenna temperature can be calculated from the noise equivalent 
power NEP and the sampling frequency Vsampiing = 200 Hz via: 

2 NEP -v/Vsampling 

= , ; {HFT) (7) 

Alternatively, for Planck's LF/-receivers, the rms-fluctuations 
(Tn in antenna temperature are given by: 

-^Av/Vsnniplmg 

Values for r„oise and NEP can be obtained from Planck's, simu- 
lation pipeline manual. The resulting effective noise level for all 
Planck channels for a single observation of a pixel is given in Ta- 
ble Q The formulae and respective parameters are taken from the- 
planck simulation manual, available via Planck's LiveLink. 

The rms-fluctuations ctn in antenna temperature have to be 
scaled by ynjet (assuming Poissonian statistics), where n^et denotes 
the number of redundant receivers per channel, because they pro- 
vide independent surveys of the microwave sky. 

From simulated scanning paths it is possible to derive an expo- 
sure map using the siininission- and multimod-utilities. An ex- 
ample of such an exposure map in the vicinity of the North ecliptic 
pole is given in Fig.Q Using the number of observations ii^bs per 
pixel, it is possible to scale down the noise amplitudes by ^fn^s 
and to obtain a realistic noise map for each channel. Here, we apply 
the simplification that all detectors of a given channel are arranged 
collinearly. In this case, the exposure maps will have sharp transi- 
tions from well-observed regions around the ecliptic poles to the 
region around the ecliptic equator. In real observations these tran- 
sitions will be smoothed out due slight displacements of the optical 
axes among each other which causes the effective exposure pattern 




Figure 1. Exposure map (side length ^ 70°) of a single v = 353 GHz- 
receiver at the North ecliptic pole in logarithmic shading: The displacement 
of the receiver with respect to the optical axis causes the observational rings 
not to overlap exactly at the pole, but gives rise to the lozenge-shaped pat- 
tern in the sky-coverage map. On average, the pixels inside the lozenge are 
observed roughly 100 times, compared to ~ 20 times outside. Pixels on the 
edges of the lozenge are observed a few thousand times. The best observed 
pixels are situated on the tips of the lozenge, where values as high as 2 ■ 10"* 
are attained. The numbers correspond to the planned mission Hfetime of 1 
year. The faint diagonal tangential line on the left side is caused by 2008's 
being a leap year: The mapping of the entire sky would be completed in 365 
days, but there is an additional day available. 



to be a superposition of rotated and distorted single-receiver expo- 
sure patterns. 



4 FOREGROUND EMISSION COMPONENTS 

The observation of the CMB and of SZ-clusters is seriously im- 
peded by various Galactic foregrounds and by the thermal emission 
of celestial bodies of our Solar system. In order to describe these 
emission compo nents, template maps from microwave surveys are 
used. Bouchet & Gisper3 il999h give a comprehensive review for 
the foreground components relevant for the Planck mission. As 
foreground components we include thermal emission from dust in 
the Galactic plane (Sect. l4m . Galactic synchrotron (Sect. 14. 21 and 
free-free emission (Sect. B3t , line emission from rotational tran- 
sitions of carbon monoxide molecules in giant molecular clouds 
(Sect. I4.4> . sub-millimetric emission from planets (Sect. I4.5> and 
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Figure 2. Sky map of the antenna temperature increase caused by dust 
emission in the v = 100 GHz-channel: The shading is proportional 
to arsinh(7'A(v = 100 GHz)/|aK). Ecliptic coordinates have been cho- 
sen. This map has b een derived from the dust-template map provided by 
ISchlegeletal1ll998t) . 



Figure 3. Sky map of the antenna temperature increase caused by syn- 
chrotron emission in the v = 100 GHz-channel in ecliptic coordinates: The 
shading is proportional to arsinh Cr^C^ = 100 GHz)/ia-K). The survey un- 
dertaken bv lHaslam et alJ il982h was used to construct this template. 



from minor bodies of the Solar system (Sect. l4.6V Foreground com- 
ponents omitted at this stage are discussed in Sect. l4.7l 

In this work, no attempt is made at modelling the interactions 
between various foreground components because of poorly known 
parameters such as the spatial arrangement along the line-of-sight 
of the emittin g and absorbing c omponents. Exemplarily, the reader 
is referred to Finkbeiner (2003), where the absorption of Galactic 
free-free emission by dust is discussed. 

4.1 Galactic dust emission 

At frequencies above ~ 100 GHz, the thermal emission from dust in 
the disk of the Milky Way is the most prominent feature in the mi- 
crowave sky. Considerable effort has been undertaken to model the 
thermal emission from Galactic dust (Schleeel etal. 19971 (1993: 
iFinkbeiner et al. ,1999- .2000) . The thermal dust emission is re- 
stricted to low Galactic latitudes and the thin disk is easily dis- 
cernible. 

The input template map (see Fig.|2} is derived from an obser- 
vation at a wavelength of A = 100 \im, i.e. vq = 3 THz. Its ampli- 
tudes Adust are given in MJy/sr, which are extrapolated to the actual 
frequency channels of Planck using a two-component model sug- 
gested by C. Baccigalupi (personal communication). Despite the 
fact that the dust is expected to spread over a large range of tem- 
peratures, the model reproduces the thermal emission remarkably 
well. This model yields for the flux SdustC"): 



S dust(v) ^ 



B(v,Ti) + f2-{^f B(v,T2) 



(9) 



f,qB(vo,Ti) + f2B(vo,T2) 

The choice of parameters used is: /i = 0.0363, /a = 1 - /i, = 
1.67, Q-2 = 2.70, q = 13.0. The two dust temperatures are Ti = 
9.4 K and T2 = 16.2 K. The function B{v, T) denotes the Planckian 
emission-law: 
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exp(hv/ksT) - 1 



(10) 



4.2 Galactic synchrotron emission 

Relativistic electrons of the interstellar medium produce syn- 
chrotron radiation by spiralling around magnetic field lines, which 
impedes CMB observations most strongly at frequencies below 
100 GHz. The synchrotron emission reaches out to high Galactic 



latitude and is an important ingredient for modelling foreground 
emission in microwave observations. An all-sky survey a t an ob 



serving frequ ency of 408 MHz has been compiled by HaslMi^^l 
^ 1981', '1982^ and adopted for usage with Planck by Giardino et al 



( 2002) (see Fig.j^. The average angular resolution of this survey is 

0;85 (FWHM). 

Recent observations with WMAP jBennett et'ai]|2003h indi- 
cate that the spectral slope of the synchrotron emission changes 
dramatically from 7 = -0.75 at frequencies below 22 GHz to 
y = -1.25 above 22 GHz. Theoretically, this may be explained by 
a momentum-dependent diffusion coefficient for cosmic ray elec- 
trons. In order to take account of this spectral steepening, the am- 
plitudes Asynchio are multiplied with a prefactor in order to obtain 
the synchrotron fluxes at v = 22 GHz. This value is then extrap- 
olated to Planck's observing frequencies with a spectral index of 
7 = -1.25: The amplitudes Asy„chro of the input map are given in 
units of MJy/sr, and for the flux 5sy„chro(''') one thus obtains: 
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(11) 



Here, the fact that the synchrotron spectral index shows signif- 
icant variations across the Milky Way due to varying magnetic field 
strength is ignored. Instead, a spatially constant spectral behaviour 
is assumed. 



4.3 Galactic free-free emission 

The Galactic ionised plasma produces free-free emission, which is 
an important source of contamination in CMB observations, as re- 
cently confirmed by Bennett et alj 1200 3h in WMAP observations. 
Aiming at modelling the free-free emission at microwave frequen - 
cies, we rely on an //^.-template provided by iFinkbeineJ i2003h . 
Modeling of the free-free emission component on the basis of an 
/f„//?/ia- template is feasible because both emission processes de- 
pend on the emission measure Jnjd/, where n^ is the number den- 
sity of electrons. This template is a composite of three //^-surveys 
and is because of its high resolution (on average 6'.0 FWHM) par- 
ticularly well suited for CMB foreground modelling. The morphol- 
ogy of the free-free map is very complex and the emission reaches 
out to intermediate Galactic latitude. 

For relating //j,-fluxes A^^ given in units of Rayleighs to the 
free-free signal's ant enna temperature Tfrec-frec measured in Kelvin, 
IValls-Gabaujjl998t) gives the formula: 
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Figure 4. Sky map of the antenna temperature increase caused by free- 
free emission in the v = 100 GHz-channel: The shading is proportional to 
arsinh(7'yi(v = 100 GHz)/ia-K). Ecliptic coordinates have been c hosen. This 
map h as been derived from the //„-template map provided bv IpinkbeineJ 
l2003l). 
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Tp denotes the plasma temperature and is set to 10"* K in this work. 
An approximation for the Gaunt factor gg valid for microwave fre- 
quencies i n the range Vn «; v <e: ksT/h (vp is the plasma frequency) 
IS given bv lFinkbeineji2003h : 
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(13) 



where e and denote electron charge and mass (in Gaussian units) 
and jE - ^.^1121 is the Euler constant. The contribution of frac- 
tionally ionised helium to the free-free emissivity as well as the 
absorption by interstellar dust has been ignored because of its be- 
ing only a small contribution in the first case and because of poorly 
known parameters in the latter case. The antenna temperature can 
be converted to the free-free flux 5frec-frec(i') by means of: 



'Sfree-free(l') — 2— • kgT . 



(14) 



Concerning the free-free emission, there might be the possi- 
bility of an additional free-free component uncorrelated with the 
/f„-emission. This hot gas, however, should emit X-ray line radia- 
tion, which has not been observed. 



4.4 CO-lines from giant molecular clouds 

In a spiral galaxy such as the Milky Way, a large fraction of the 
interstellar medium is composed of molecular hydrogen, that re- 
sides in giant molecular clouds (GMC), objects with masses of 
lO'' - 10' Mq and sizes of 50 - 200 pc. Apart from molecular 
hydrogen, the GMCs contain carbon monoxide (CO) molecules 
in significant abundance. The rotational transitions of the CO 
molecule at 115 GHz and higher harmonics thereof constitute a 
source of contamination for all P/afic/r/fFZ-channels. An extensive 
search for ato mic and molecular transition lines was undertaken by 
iBennettet ai] fl994^ with the FIRAS instrument onboard COBE. 

The CO-contamination is inodelled by employing a mosaic of 
CO-surveys assembled by Dame et al. 1 1996, 2001). It shows the 
velocity-integrated intensity of the transition from the first excited 
state (7 = 1) to the ground state {] = 0) close to the Galactic plane 
{jj < 5°), and additionally comprises a few CO clouds at higher 
Galactic latitude, as well as the Large Magellanic Cloud and the 



Figure 5. Sky map of the increment in antenna temperature due to CO-line 
emission in the v = 100 GHz-channel in ecliptic coordinates: The shading 
is proportional to arsinh(rA(v = 100 GHz)/|aK). The maps shows the ro- 
tational transition of the CO molecule from the first excited s tate in to the 
ground state at voni = 115 GHz as derived bv lPameetal] 1200 ll) for a 
temperature of Tqq = 20 K. 



Andromeda galaxy M 31. Due to the composition of the map, the 
angular resolution is not uniform, but the best resolution of ^ T.5 
is reached for a large area around the Galactic plane. 

From this map, it is possible to derive the line intensities of the 
higher harmonics, assuming thermal equilibrium: The frequency 
V for a transition from a state of rotational quantum number / to 
a state with quantum number 7 -I- 1 of the CO molecule follows 
from elementary quantum mechanics: The rotational energy of a 
molecule with moment of inertia 6 and angular momentum / is 
Emt = /20 = fr ■ J( J +1) 126. In the last step the quantum number 
/ was introduced. For the transition energy between two subsequent 
rotation levels, one obtains: 



vj„j+i = 2ec ■ (y + 1) = 1 15 GHz ■ (7 + 1), 



(15) 



where Q = hj^n^cO is a measure of the inverse moment of inertia 
of the molecule and c denotes the speed of light. Thus, the spectrum 
consists of equidistant lines. The relative intensities of those lines 
is given by the ratio of their occupation numbers ^y: 



, Qhc 

Xj = {2J +\)-&x^\—^—J{J+\) 



(16) 



i.e. the relative line intensities of two consecutive lines is 

given by: 



Xj 



2J + 3 
2J+ 1 



• exp 



2Qhc 
kgT CO 



(J+l) 



(17) 



Xj IS, detemined by a statistical weight (27 +1) reflecting the de- 
generacy of angular momentum and a Boltzmann factor. For the 
determination of line intensities thermal equilibrium is assumed, 
common estimates for the temperature inside GMCs are Tqo = 10- 
30 K. For the purpose of this work, we choose Tco = 20 K. From 
the brightness temperature Ta one obtains the CO-flux 5co-ime(v) 
by means of the following equation: 



5C0-Iine(v) = 2— • ksTAi^) ■ p(v - Vj^j+i), 



(18) 



where the line shape p(v - vj^j+i ) is assumed to be small in com- 
parison to Planck's frequency response windows such that its actual 
shape (for instance, a Voigt-profile) is irrelevant. 
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4.5 Planetary submillimetric emission 

Planets produce infra-red and sub-millimetric radiation by absorb- 
ing sunlight and by re-emitting this thermal load imposed by the 
sun. The investigation of the thermal properties of Mars, Jupiter and 
Saturn has be en the target of several space missions ( Goldin et alj 
1 19971 iGriffin e t al. 1986, to name but a few). For the descrip- 
tion of the submillimetric thermal emission propertie s of plan- 
ets, an extension to the Wright & Odenwald model ( Wri ghJl976t 
iNeugebauer et alll97ll) was used. The orbital motion of the plan- 
ets is sufficiently fast such that their movements including their 
epicyclic motion relative to the Lagrangian point Li, Planck's ob- 
serving position, has to be taken into account. All planets are im- 
aged twice in approximate half-year intervals due to Planck's scan- 
ning strategy, while showing tiny displacements from the ecliptic 
plane because of the Lissajous-orbit of Planck around Li and their 
orbital inclinations. 

The heat balance equation for a planet or asteroid reads as: 



E + F +W = P^, 



P. 



Libsoiption 



= I + R, 



(19) 



where E denotes the heat loss by thermal emission (i.e. the signal 
for Planck), F the heat flux outward from the interior of the planet, 
W is the heat lost by conduction to the planet's atmosphere, / is the 
Solar radiation absorbed and R is the heating of the planet caused 
by the back- scattering of radiation emanating from the surface of 
the planet by the atmosphere. The definition of these quantities is 
given by eqns. <20> through <24> : 
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(20) 
(21) 
(22) 

(23) 
(24) 



Here, e is the surface emissivity of the planet, cr is the Stefan- 
Boltzmann constant, Tpianct is the planet's temperature, k the co- 
efficient of heat conduction, A the planet's bolometric albedo, G 
the Solar constant (i.e. the energy flux density of Solar irradiation 
at the earth's mean distance), r the distance of the planet to the sun 
in astronomical units, t the planet's rotation period and 6' the ge- 
ographical latitude of the radiation absorbing surface element. The 
temperature distribution in the interior of the planet at radial posi- 
tion X is controlled by the heat conduction equation: 



97-, 



planet 



planet 



(25) 



9? 9x2 ' 

with the specific heat per unit volume c. 

In our model, the heat loss R of the planet's surface due to con- 
duction to the planet's atmosphere is taken to be a constant fraction 
of the heat flux F outward from the interior of the planet, the con- 
stant of proportionality being k, for which we assumed k — 0.1. 
Similarly, the heat gain by back-scattering radiation by the atmo- 
sphere R was assumed to be a constant fraction y of the local noon 
Solar flux /„,ax, where y was taken to be y = 0.01. The system of 
differential eqns. j20> - J24> dependent on time t and on Solar dis- 
tance r constitutes a heat conduction problem with periodic excita- 
tion (by the planet's rotation). Thus, the heat balance of the planets 
is modelled by periodic solutions of the Laplacian heat conduction 




Figure 6. Sky map of increment in antenna temperature Ta produced 
by planets and asteroids in the v = 30 GHz channel at 33'4 reso- 
lution (FWHM). The colour coding is proportional to arsinh(7'^(v = 
100 GHz)/nK). The asteroids reach ecliptic latitudes of \I3\ i, 30°. The pro- 
nounced signals are produced by planets, which (due to Planck's scanning 
strategy) appear twice. The observable planets comprise (from left to right) 
Saturn, Mars, Uranus, Neptune and Jupiter. The epicyclic motion of Mars 
is sufficiently fast to counteract the parallactic displacement such that it ap- 
pears only once. 



differential equations. It was solved iteratively by applying Laplace 
transforms with periodic boundary conditions. The integration over 
the planet's surface then yields the radiation flux. In the calculation, 
we addressed rocky and gaseous planets differently with respect to 
their thermal properties. Furthermore, the giant gaseous planets are 
known to have internal sources of heat generation, which also has 
been taken account of. 

The brightest point source in the microwave sky due to the 
planetary thermal emission is Jupiter, causing an increase in an- 
tenna temperature of Tjupiie, = 93.6 mK in the v = 100 GHz- 
channel, followed by Saturn with Tsatum = 15.0 mK. All outer plan- 
ets apart from Pluto will be visible for Planck. Estimates show that 
even Galilean satellites Ganymede, Callisto, lo and Europa and Sat- 
urn's moon Titan are above the detection threshold of Planck, but 
they are outshone by the stray-light from Jupiter and Saturn, re- 
spectively and for that reason not included in our analysis. 

Due to the planet's being point sources, their fast movement 
and their diverse surface temperatures it is not feasible to produce a 
template and extrapolate the fluxes with a common emission law to 
/'/ancA;-frequencies. Instead, flux maps have been produced directly 
for each of the nine Planck-channels separately taking account of 
the planetary motion, the solution of the heat balance equation laid 
down above and the finite beam-width. The analogous holds for 
asteroids, that are covered by the next chapter. 



4.6 Submillimetric emission from asteroids 

Asteroids and minor bodies of the Solar system are easily observed 
by infrared satellites such as ISO and possibly by sub-millimetric 
observa tories (Mii ller 2001; Midler & Laaerros 1999). An estima- 
tion bv lCremonese et alj ^2002) shows that a large number of as- 
teroids (~ 400) should yield signals detectable by Planck. The or- 
bital motion of all asteroids is fast enough to cause double detec- 
tions at different positions in the sky separated by half a year due 
to Planck's scanning strategy. In contrast to planets, asteroids are 
not well restricted to the ecliptic plane and appear up to ecliptic 
latitudes of/3i 30°. 

The thermal emission properties of asteroids are well un- 
derstood (for a comprehensive and detailed review, see iLagerrosI 
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Il996albl 1 19971 Il998h such that asteroids have been used for cal- 
ibrating det ectors (e.g. the ISO mission, c.f. iMiiller & Lagerro^ 
1199811200 2^^ and for determining beam shapes. The thermal model 
used for describing the submillimetric emission by asteroids is the 
same extension of the Wright & Odenwald model as for rocky plan- 
ets. However, additional features that had to be incorporated was 
the beamed emission due to surface roughness. Furthermore, in the 
system of differential eqns. <20> - <24> terms W and R were ne- 
glected due to the absence of atmospheres in asteroids. 

Information about the diameter and albedo was derived us- 
ing the HG-magnitude system in case of asteroids for which those 
quanti ties are unkn own, otherwise literature values were taken 
(from lMoor'3l200l and lAU's Minor Planet Centre ^). For the 
description of the rotation period, an empirical relation that ex- 
presses the rotation period as a function of mass was used in the 
cases where the rotation period is unknown. The brightest sources 
include Ceres (Tccics = 19.7 |j.K), Pallas (Tf jiaa = 7.2 |i.K), Vesta 
(T'vesta = 6.7 ii.K) and Davida (rD^vida = 2. 1 ii.K). The temperatures 
stated are antenna temperatures measured in the v = 100 GHz- 
channel at the brightness maximum. 

Our simulation shows that the number of detectable asteroids 
is overestimated by Cremonese et al. (2002), who did not take the 
expected observation geometry and detector response into account. 
Typical surface temperatures of asteroids are of the order of 150 K, 
and therefore, Planck is observing their thermal emission in the 
Rayleigh-Jeans regime. For that reason, the number of detectable 
asteroids increases with observing frequency. For our sample of 
5 • lO'* asteroids of the Minor Planet Centre's catalogue, we find 
a couple of asteroids at v = 30 GHz, a few tens of asteroids at 
V = 100 GHz and up to 100 asteroids in the highest frequency band 
at V = 857 GHz. Approximately 1200 asteroids will have fluxes 
above half of Planck's single-band detection limit estimated for 
ideal observation conditions and thus they constitute an abundant 
population of point sources that possibly hampers the detection of 
SZ-clusters. 

The prediction of comets is very uncertain for the years 2007 
through 2009: Many comets are not detected yet, non-active comets 
are too faint with few exceptions and the coma thermal emission 
features of active comets is very complex. For these reasons, they 
have been excluded from the analysis. 

4.7 Future work concerning Planck\ foregrounds 

Foreground components not considered so far include microwave 
point sources, such as infra-red galaxies and microwave emitting 
AGNs. The emission of infra-red galaxies is associated with ab- 
sorption of star light by dust and re-emission at longer wavelengths. 
Galaxies with ongoing star formation can have large fractions 
(~ 90%) of their total emission at infra-red wavelengths, compared 
to about one third in the case of local galaxies. The integrated emis- 
sion from unresolved infra-r ed galaxies accounts for the cosmic 
infra- red background (GIB) ( Puget et alJll99d: iLagache & Puged 
l2OO0h . the fluctuations of which ar e impeding SZ-ob servations at 
freq uencies above v ^ 100 GHz ( Aghanin^^lJ|200^. 

iLagachd j2003) and iWhite & Maiu mdar' flOOf) have esti- 
mated the number counts of unresolved infra-red jjalaxies at 
Planck frequencies, which was used by Aghanim et alJ j2004) in 
order to estimate the level of fluctuation in the Planck-beam. In the 
easiest case, the sources are uncorrelated and the fluctuations obey 

^ Ihttp ://cfa-www. harvard ■ edu/cfa/ps/mpc.html| 



Poissonian statistics, but the inclusion of corre lations is expecte d 
to boost the fluctuations by a facto r of ~ 1.7 iSong et aljEooS) . 
According to Aghani m et alJi2004) . the resulting fluctuations vary 
between a few 10^ Jy/sr and 10^ Jy/sr, depending on observing 
channel. A proper modeling would involve a biasing scheme for 
populating halos, the knowledge of the star formation history and 
template spectra in order to determine the K-corrections. 

AGNs are another extragalactic source of submillimetric emis- 
sion. Here, sychrotron emission is the radiation generating mech- 
anism. The spectra show a variety of functional behaviours, with 
spectral indices a generally ranging from -1 to -0.5, but sources 
with inverted spectra a > are commonplace. This variety makes it 
difficult to extrapolate fluxes to observ ing frequencies of CMB ex- 
periments. Two studies jToffolatti et"aU1998; ,Sokasian et al..20oil) 
have estimated the fluctuations generated by radio emitting AGNs 
at SZ-frequencies and found them to amount to 10^ - lO"* Jy/sr. 
However, AGNs are known to reside in high-density environments 
and the proper modelling would involve a (poorly known) biasing 
scheme in order to assign AGN to the dark matter halos. Apart from 
that, one would have to assume spectral properties from a wide 
range of spectral indices and AGN activity duty cycles. Therefore, 
the study of extragalactic sources has been omitted from this anal- 
ysis. 

Yet another source of microwave emission in the Solar sys - 
tem is the zodiacal light jLeinert et ai]|2002l: iReach et alj|2003 j) . 

Modelling of this emission component is very difficult due to the 
Lissajous-orbit of Planck around the Lagrangian point L2. The 
disk of interplanetary dust is viewed under varying angles depend- 
ing on the orbital period and the integration over the spatially 
non-uniform emission features is very complicated. Reach et al] 
(2003b) have investigated the thermal emission by interplanetary 
dust from measurements by ISO and have found dust temperatures 
of rjodiiicai = 250 - 300 K and fluxes on the level of ^ 10^^ Jy/sr, 
i.e. the equilibrium temperature is separated by two orders of mag- 
nitude from the CMB temperature, which means that the intensi- 
ties are suppressed by a factor of ~ 10"* due to the Rayleigh-Jeans 
regime of the zodiacal emission in which Planck is observing and 
by a factor of 10' due to Planck's narrow beams. From this it is con- 
cluded that the emission from zodiacal light is unlikely to exceed 
values of a few ~ (j-Jy in observations by Planck which compares 
to the fluxes generated by faint asteroids. Thus, the zodiacal light 
constitutes only a weak foreground emission component at submil- 
limetric wavelengths and can safely be neglected. 



5 SIMULATING SZ-OBSERVATIONS BY Planck 

The simulation for assessing Planck's SZ-capabilities proceeds in 
four steps. Firstly, all-sky maps of the thermal and kinetic SZ- 
effects are prepared, the details of map-construction are given in 
Sect. 15.11 Secondly, a realisation of the CMB was prepared for 
the assumed cosmological model (Sect. l?!2t . The amplitudes were 
co-added with the Galactic and ecliptic foregrounds introduced 
in the previous section, subsequently degraded in resolution with 
Planck's beams (Sect. 15. 3> . Finally, uncorrelated pixel noise as 
well as the emission maps comprising planets and asteroids were 
added. In the last section, cross-correlation properties of the vari- 
ous astrophysical and instrumental noise components are discussed 
(Sect.|53J. 

At this stage it should be emphasised that we work exclu- 
sively with spherical harmonics expansion coefficients Of,,, of the 
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flux maps. The expansion of a function a(0) into spherical harmon- 
ics Y"'{ff) and the corresponding inversion is given by: 



dfi a{e) ■ Y'I'iO)' and a{e) = J]J] «f,„ • Y^(9). (26) 



Here, df2 denotes the differential solid angle element. For reasons 
of computational feasibility, we assume isotropic spectral proper- 
ties of each emission component, i.e. the template map is only pro- 
viding the amplitude of the respective emission component, but the 
spectral dependences are assumed to remain the same throughout 
the sky. While this is an excellent approximation for the CMB and 
the SZ-effects (in the non-relativistic limit), it is a serious limitation 
for Galactic foregrounds, where e.g. the synchrotron spectral index 
or the dust temperatures show significant spatial variations. 

Adopting this approximation, the steps in constructing spher- 
ical harmonics expansion coefficients (S [„,)yg of the flux maps 
S (0, v) for all Planck channels consist of deriving the expansion 
coefficients of the template, converting the template amplitudes to 
flux units, extrapolate the fluxes with a known or assumed spectral 
emission law to Planck's observing frequencies, to finally convolve 
the emission law with Planck's frequency response window for 
computing the spherical harmonics expansion coefficients of the 
average measured ffux (5f„,)v„ at nominal frequency vo by using 
eqn. 1271 . 



(S Im) 



/dv5 



fm(v)/?^„(v) ^ 1^ 



/dvi?.„(v) 



kRTr, 



(27) 



Here, 5f,„(v) describes the spectral dependence of the emis- 
sion component considered, and /?vo(v) the frequency response of 
Planck's receivers centered on the fiducial frequency vg. Assum- 
ing spatial homogeneity of the spectral behaviour of each emis- 
sion component it is possible to decompose 5f,„(v) into Stm(v) = 
^(v) • Of,,,, i.e. a frequency dependent function q(v) and the spherical 
harmonics expansion coefficients a^,,, of the template describing the 
morphology. This is possible due to the fact that the decomposition 
eqn. <26> is linear. Additionally, eqn. <27> gives the conversion from 
the averaged flux (S [,„)v in a Planck-chaxmA to antenna tempera- 
ture Tcm- 

Planck's frequency response function Rygiy) is well approxi- 
mated by a top-hat function: 



1, V 6 [vo - Av, Vf) + Av] 
0, V ^ [vo - Av, Vo + Av] 



(28) 



The centre frequencies vq and frequency windows Av for Planck's 
receivers are summarised in Table. Q In this way it is possible to 
derive a channel-dependent prefactor relating the flux expansion 
coefficients <5f„,),,„ to the template expansion coefficients The 
superposition of the various emission components in spherical har- 
monics and the determination of response-folded fluxes is most 
conveniently done using the almmixer-utility of Planck's simu- 
lation package. 

5.1 SZ-map preparation 

For constructing an all-sky Sunyaev-Zerdovich map, a hybrid ap- 
proach has been pursued. Due to the SZ-clusters being detectable 
out to very large redshifts, due to their clustering properties on very 
large angular scales, and due to the requirement of reducing cos- 
mic variance when simulating all-sky observations as will be per- 
formed by Planck, there is the need for very large simulation boxes, 
encompassing redshifts of z ^ 1 which corresponds to comoving 
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ecliptic longitude .1 [deg] 

Figure 7. Detail of the thermal Comptonisation map: A 2° X 2° wide cut- 
out centered on the ecliptic coordinates (/i,/3) = (45°, 45°) is shown. The 
smoothing imposed was a Gaussian kernel with A6 = 2^0 (FWHM). The 
shading indicates the value of the thermal Comptonisation y, which is pro- 
portional to arsinh(10* ■ y). This map resulted from a projection on a Carte- 
sian grid with mesh size ~ 14", i.e. no HEALPix pixelisation can be seen. 
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Figure 8. Detail of the kinetic Comptonisation map: A 2° X 2° wide cut- 
out centered on the same position as Fig.Q i.e. at the ecliptic coordinates 
(/!,/?) = (45°, 45°) is shown. The smoothing imposed was a Gaussian kernel 
with AS = 2^0 (FWHM). The kinetic Comptonisation w is indicated by the 
shading, being proportional to arsinh(10' ■ w). 



scales exceeding 2 Gpc. Unfortunately, a simulation incorporating 
dark matter and gas dynamics that covers cosmological scales of 
that size down to cluster scales and possibly resolving cluster sub- 
structure is beyond computational feasibility. For that reason, two 
simul ations have been combined: The H ubble-volume simulation 
I Jenk ins et all200ll:IColberg et alj200ol) . and a smaller scale simu- 
lation including (adiabatic) gas physics by White et al. ( 2002) per- 
forme d with GADGET I Springel et al. 2001; Springel & Herna ui's3 
l2002t) . 

All-sky maps of the SZ-sky were constructed by using the 
light-cone output of the Hubble-volume simulation as a cluster cat- 
alogue and template clusters from the small-scale gas-dynamical 
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dimensionless frequency x = hv/iksTcMB) 



Figure 9. Frequency dependence of the thermal (thick hnes) and of the ki- 
netic SZ-flux (thin Hnes), for ideal i5-like frequency responses (solid lines), 
for a top-hat window function with a relative width of 10% corresponding to 
Planck's LF/-instrument and for a top-hat window function with a relative 
width of 16.7%, as planned for Planck's //^/-instrument. The fluxes stated 
correspond to the integrated Comptonisation of J/ = "W = 1 arcmin" . The 
vertical lines indicate the centre frequencies of Planck's receivers. 



simulation. In this way, the sky-maps contain all clusters above 
5 • lO'^M(T)/h out to redshift z = 1.48. The analysis undertaken 
by^artelmann (2001) gives expected mass and redshift ranges for 
detectable thermal SZ-clusters, which are covered completely by 
the all-sky SZ-map presented here. The maps show the correct 2- 
point halo correlation function, incorporate the evolution of the 
mass function and the correct distribution of angular sizes. 

Furthermore, they exhibit cluster substructure and deviations 
from the ideal cluster scaling relations induced by the departure 
from spherical symmetry. The velocities used for computing the 
kinetic SZ-effect correspond to the ambient density field. The map 
construction proces s and the properties of the resulting map are in 
detail described in iSchafer et alj i2004) . Visual impressions of the 
SZ-maps are given by Figs.0and|8l 

The fluxes generated by the thermal SZ-efFect Sj/(x) and of 
the kinetic SZ-efFect S-\y(x) are given by eqns. and QUI , 
respectively. The dimensionless frequency is defined as x = 
hv/ikgTcMB) and the flux density of the CMB is given by 5o = 
= 22.9 Jy/arcmin^: 



Si/(x) 
S-wix) 



• exp(j) 



(exp(x) - 1)2 



exp(x) + 1 

c 

exp(x) - 1 



■ exp(x) 



(exp(x)- 1)2 



(29) 
(30) 



TableQ summarises the fluxes 5 j/ and S-w and the correspond- 
ing changes in antenna temperature Ty and for the respective 
Comptonisation of J/ = = 1 arcmin^ for all /'/ancA:-channels. 

Fig.|9|shows how the frequency dependence of the SZ-signal 
is altered by Planck's relatively broad frequency response func- 
tions. The relative deviations of curves in which the frequency win- 
dow has been taken into account to the unaltered curve amounts to 
5 . . . 15%, depending on observation frequency. 



5.2 CMB-map generation 

The angular power spectrum Cf is comput ed for a flat ACDM- 
cosmology using the CMBfast code by ISeliak & Zaldarriagal 
(■1996*). In addition to the cosmological parameters being already 
given in Sect. Q we use adiab atic initial co nditions, set the CMB 
monopole to Tcmb = 2.725 K ^Mather et al . 1999) and the primor- 
dial He-mass fraction to X^f. = 0.24. The reionisation optical depth 
T was set to t = 0.17 and the reio nisation redshift was taken to 
be Ziaon = 20 teennett et 31.11200 3^. The angular power spectrum 
of the CMB is normalised to COBE data. With the spectrum of 
Cf-coeflicients, a set of a^,,, -coefficients was synthesi s ed by using 
the synalm code based on synfast bv iHivon et alj il99^ . The 
factors for converting the af„i -coefficients of the CMB map show- 
ing the thermodynamic temperature and to the corresponding fluxes 
for each channel were then derived by convolution of the Planckian 
emission law eqn. <3H . 



S cmb(v) - Sq ■ 



(31) 



exp(j) - 1 ' 

with Planck's frequency response function eqns. 127 > and <28> . 
Again, 5o = 22.9 Jy/arcmin^ is the energy flux density of the CMB. 



5.3 Preparation of simulation data sets 

The expansion coefficients of the flux maps are multiplied with the 
respective beam's fcfo-coefficients in order to describe the finite an- 
gular resolution. After that, expansion coefficients of the pixel noise 
maps and those of the planetary maps have been added. In total, 
three atlases consisting of nine flux (5f,„)^„-sets belonging to each 
of Planck's channels with fiducial frequency vo have been com- 
piled: 

• The reference data set is a combination of the CMB, the SZ- 
maps and the instrumental noise maps. They should provide the 
cleanest detection of clusters and the measurement of their proper- 
ties. Apart from the inevitable instrumental noise, this data set only 
contains cosmological components. In the remainder of the paper, 
this data set will be refered to as COS. 

• The second data set adds Galactic foregrounds to the CMB, 
the SZ-maps and the instrumental noise map. Here, we try to as- 
sess the extend to which Galactic foregrounds impede the SZ- 
observations. Thus, this data set will be denoted GAL. 

• In the third data set the emission from bodies inside the So- 
lar system was included to the CMB, the SZ-maps, the Galactic 
foregrounds and the instrumental noise. Because of the planets and 
asteroids being loosely constrained to the ecliptic plane, this data 
set will be called ECL. 

An example of a synthesised map showing the combined 
emission of the SZ-clusters and all Galactic and ecliptic compo- 
nents including neither CMB fluctuations nor instrumental noise 
at a location close to the Galactic plane is given by Fig. 1101 The 
observing frequency has been chosen to be v = 143 GHz, corre- 
spondingly, the map has been smoothed with a (Gaussian) beam of 
A0 = 7:1 (FWHM). 



5.4 Planck-channel correlation properties 

In this section the auto- as well as the cross-correlation properties 
of the various foregrounds in different Planck-channels are studied. 
The cross power specta, defined formally by eqn. iili are deter- 
mined by using: 
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Figure 11. Power spectra in various P/a«c/;-channels: v = 30 GHz (solid), 
V = 143 GHz (dashed), y = 353 GHz (dasli-dotted) and v = 857 GHz 
(dotted) for COS data set (tliin line), the GAL data set (medium line) and the 
ECL data set (thick line). 



Figure 10. A 50° x 50° wide composite centered on the ecliptic coordi- 
nates (/!,j8) = (-85°, 0°), i.e. close to the Galactic centre for Planck's 
V = 143 GHz-channel. The shading is proportional to aismh(TA(y = 
143 GHz)/mK). The map is smoothed with the corresponding beam of di- 
ameter A8 = T.l (FWHM). SZ-clusters are observed in absorption in this 
channel and are discernible by eye even at close proximity (b 20°) to the 
Galactic plane. For clarity, the CMB fluctuations as well as the instrumental 
noise have been excluded. The two point sources on the ecliptic equator are 
twin detections of Jupiter. 



2e + 



■(S(,„r„ 



(32) 



From this definition, the auto-correlation spectra are obtained by 
setting vi = V2, i.e. Cr,v = Cr.w The band-pass averaged fluxes 
(S f„,)y are defined in eqn. <27t . In Fig. 1111 the power spectra are 
shown for the v = 30 GHz-, v = 143 GHz-, v = 353 GHz- and the 

V = 847 GHz-channels. The spectra have been derived including 
various Galactic and ecliptic noise components in order to study 
their relative influences. For visualisation purposes, the spectra are 
smoothed with a moving average filter with a filter window com- 
prising 11 bins. 

Distinct acoustic peaks of the CMB are clearly visible in the 
clean COS data sets, but are overwhelmed by the Galactic noise 
components. At small scales, i.e. high multipole order £, differ- 
ences between the GAL and ECL data sets become apparent, the 
latter showing a higher amplitude. The (single) acoustic peak mea- 
surable in the v = 33 GHz channel is shifted to larger angular scales 
due to the coarse angular resolution of that particular channel. The 

V = 857 GHz-curve of the COS data set behaves like a power law 
due to the fact that the CMB is observed in the Wien-regime and 
is consequently strongly suppressed, such that the angular power 
spectrum is dominated by uncorrelated pixel noise. 

Fig. 1121 shows exemplarily a couple of cross power spectra. 
The cross-correlation spectra derived for the COS data set nicely 
shows the CMB power spectrum if two neighboring channels close 
to the CMB maximum are chosen, but the correlation is lost in two 
widely separated channels. This is especially the case if one con- 
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Figure 12. Cross-correlations: The spectra Cf_vj=i43 ghz,v2=545 GHz (solid 
line) and Cf,vi=l43 GHz,v2=2l7 GHz (dashed line) are contrasted for the 
COS data set (thin lines) versus the GAL data set (thick line). Further- 
more, the spectrum Cf,yj=545 gHz,v-)=857 GHz (dash-dotted line) as well as 
Cf,v|=3o GHz,v2=44 GHz (dotted line) is shown as derived from the ECL data 
set. 



siders the two lowest LF/-channels at angular scales which the re- 
ceivers are not able to resolve. In this regime the pixel noise is still 
very small and the cross-correlation spectrum drops to very small 
values. 

In order to illustrate the complexity of spectral and morpho- 
logical behaviour of the power spectra, they are given as contour 
plots depending on both the observing frequency v and the multi- 
pole order £. Fi g. 1 1 3 1 and 1 1 41 contrast the auto-correlation proper- 
ties of the different data sets. The COS data set, shown in Fig. 1131 
containing nothing but the CMB and instrumental noise apart from 
the SZ-contribution, shows clearly the acoustic oscillations with 
the first peak ai £ ^ 200 and the consecutive higher harmonics. 
They are most pronouced in the y = 100 GHz- and v = 143 GHz- 
channels. At higher multipole moments, the power spectra are dom- 
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Figure 13. Auto-correlations: The power spectral Cf y-coefficients are 
shown as a function of observing frequency v and multipole order C in the 
usual representation [{£ + l)Cf^y/2jT. The amplitudes are given in loK^ and 
the contours are linearly spaced. Note the logarithmic scaling of the fre- 
quency axis. In the data set displayed, the CMB, both SZ-effects and in- 
strumental noise are included. The first three acoustic oscillation peaks are 
clearly visible. 
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Figure 14. Auto-con'elations: The power spectral Cf_,,-coelficients are 
shown as a function of observing frequency v and multipole order f in the 
usual representation ((t + \)Ci^yl2n. The ampHtudes are given in log(|jK-) 
with logarithmically equidistant contours. In the data set displayed, the 
CMB, both SZ-elfects, all Galactic foregrounds and instrumental noise are 
included. 



inated by instrumental noise which leads to a rapid (power law) 
incline. 

Adding Galactic foregrounds yields the spectra depicted in 
Fig. ll4l Inclusion of Galactic foregrounds significantly complicates 
the picture and masks off the primary anisotropies. The spectra are 
dominated by large-scale emission structures of the Milky Way, 
most notably the emission from thermal dust that causes the spectra 
to increase with increasing frequency v. 



6 CLUSTER DETECTION BY USING 

MULTI-FREQUENCY OPTIMISED FILTERING 

One challenge in the analysis of two-dimensional all-sky surveys 
is the extraction of sources of interest which are superposed on a 
background of noise of varying morphology and spectral behaviour. 
In the presence of small-scale noise the conventional method to 
extract sources is low-pass filtering (e.g. with a Gaussian kernel) 
while wavelet analysis is most suitably applied if large scale noise 
fluctuations dominate. These methods, however, fail if the charac- 
teristic scale of the background fluctuations is comparable with the 
scale of the signal structures. Other methods have been proposed 
in order to separate difl^erent components in multifrequency CMB 
observations: They include Wiener filtering ( Tegmark & Efstathio^ 
Il99d: iBouc het & Gispert'l999: Bouchet et al. 1999), maximum- 
entropy rnethods_^Iobsoi^tal^^9^^999), Mexican-hat wavelet 
analysis jVielvaetal|l200ll:ICavon e^JjipOO), fast independent 
component analys"is~fMainoe^L'^200^, matched filter analysis 
(Tesmark & de Oliveira-Costa 1998), adaptive filtering techniques 
dSanz et alj 1200 ll : |Herranz et al. 2002b a), and non-parametric 
Bayesian approaches tPiego et alj2002.') . 

However, a comparison between these methods is difficult be- 
cause all of them assume different priors about the spatial proper- 
ties and frequency dependence. Using prior knowledge about the 
frequency dependence and statistical properties of several images 
at different frequency channels, the maximum-entropy method and 
Wiener filtering are able to separate the components of interest. 
Contrarily, wavelet analysis is well suited in order to detect com- 
pact sources. A combination of these different technique s improves 
the quality of component separation iVielva et all200ll) . Although 
component separation methods which assume a prior knowledge 
about the data are quite powerful, they yield biased or even wrong 
results in the case of incorrect or idealised assumptions about the 
data. Any error in the separation of one component propagates to 
the separation of the other components owing to normalisation con- 
straints. In particular, this is the case in non-centrally symmetric 
source profiles, oversimplified spectral extrapolations of Galactic 
emission surveys into other wavebands, variations of the assumed 
frequency dependence, or non-Gaussian noise properties the statis- 
tics of which can not fully be characterised by power spectra. Thus, 
the application of a specific component separation method is a 
trade-off between robustness and effectiveness with regard to the 
particular problem. 

Filtering techniques relying on Mexican-hat wavelets and on 
matched and scale-adaptive filters are single component separation 
methods. They all project either spatial structure or frequency prop- 
erties (within a given functional family) of the component of inter- 
est in the presence of other components acting as background in this 
context. While Mexican-hat wavelet analysis assumes Gaussian 
profiles superimposed on large scale variations of the background 
noise, the matched and scale-adaptive filter generalises to arbitrary 
source profiles and noise properties which are assumed t o be lo- 
cally homogeneous and isotropic (Sanz et al. 2001; Herra rizet alj 
E002iihft . 

This section generalises the matched and scale-adaptive fil- 
ter techniques to global spherical topologies which find ap- 
plication in all-sky surveys such as the case of Planck'?, mi- 
crowave/submillimetric survey. In addition, optimised filters for the 
detection of compact sources in single frequency all-sky observa- 
tions are derived in the appendix in a more detailed fashion. The 
proposed method aims at simultaneously localising SZ-clusters and 
measuring both their amplitudes and angular extent. It can also be 
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applied for localising microwave point sources and estimating their 
spectral properties. 

We choose the spherical filtering approach rather than tiling 
the sky with a set of two-dimensional flat maps for the following 
reasons: On the sphere, we do not have to worry about spurious 
detections as well as double detections due to overlaps in the tes- 
selation. Secondly, our approach provides a physical interpretation 
of our filter shapes in harmonic space even for the smallest multi- 
pole moments in contrast to the case of a flat map where the small- 
est wavenumbers are determined by the map size. Finally, our ap- 
proach circumvents projection failures of the noise properties such 
as stretching effects in the case of conformal mapping which would 
introduce artifical non-Gaussianity in our maps and distort profile 
shapes close to the map boundaries. 

We pursue the concept of the multi-frequency appro ach rather 



ncy appro acn r 

than the combination method (c.f. 'Herra nz et al]l2002ah . In other 
words, we filter each channel separately while taking into account 
the different cross-correlations between the different channels and 
the frequency dependence of the signal when constructing the opti- 
mised filters. This method seems to be superior to the combination 
method which tries to find a optimised combination of the different 
channels with regard to the signal-to-noise ratio of the sources and 
successively applies filters to the combined map. 

The concept is introduced and central definitions are laid down 
in Sect. 16.1 1 The concept of constructing filter kernels is outlined 
in Sect. 16. 2I Subsequently, the matched and scale-adaptive filters 
are derived for expansions of spherical data sets into spherical har- 
monics in Sect. l6.3i and Sect. 16.41 Then, the numbers of merit are 
defined in Sect. 16.51 Caveats in the numerical derivation are listed 
in Sect. 16.61 A discussion of filter kernel shapes in Sect. 16.71 for 
actual simulation data. The application of the filter kernels to our 
simulated sky maps and the extraction of the SZ-cluster signal is 
described in Sect. 16. 81 



6.1 Assumptions and definitions 

When constructing the particular filters, we assume centrally sym- 
metric profiles of the sources to be detected. This approximation is 
justified for most of the clusters of Planck's sample whose angular 
extent will be comparable in size to Planck's beams, i.e. the in- 
strumental beam renders them azimuthally symmetric irrespective 
of their intrinsic shape. Azimuthal symmetry is no general require- 
ment for the filters which can be generalised to detect e.g. elliptic 
clusters using expansions into vector rather than scalar spherical 
harmonics. 

We furthermore assume the background to be statistically 
homogeneous and isotropic, i.e. a complete characterisation can 
be given in terms of the power spectrum. This assumption obvi- 
ously fails for non-Gaussian emission features of the Galaxy or of 
the exposure- weighted instrumental noise on large angular scales. 
However, the spherical harmonics expansion of any expected com- 
pact source profile, which we aim to separate, peaks at high values 
of the multipole moment due to the smallness of the clusters where 
the non-Gaussian influence is negligible. Thus, we only have to re- 
quire homogeneity and isotropy of the background on small scales. 

In order to construct our filters, we consider a set of all-sky 
maps of the detected scalar field s^iO) for the different frequency 
channels 

sAO) = fyyy(\0 - 0o\) + ny(0), v^l,...,N, (33) 
where 6 = (&, Lp) denotes a two-dimensional vector on the sphere. 



00 is the source location, and N is the number of frequencies (re- 
spectively, the number of maps). The first term on the right-hand 
side represents the amplitude of the signal caused by the thermal 
and kinetic SZ-effect, y{\6 - 9o\) and w(\6 - 6o\), respectively, while 
the second term corresponds to the generalised noise which is com- 
posed of CMB radiation, all Galactic and ecliptic emission com- 
ponents, and additional instrumental noise. The frequency depen- 
dence of the SZ-effect is described by /,, in terms of average flux. 



fv = {S}))y and/^ = <5tv),, 



(34) 



where {S )y denotes the flux weighted by the frequency response at 
the fiducial frequency v (c.f. eqn. <27» and Sy and S'w denote the 
SZ-fluxes given by eqns. and CT . 

We expect a multitude of clusters to be present in our all-sky 
maps. In order to sketch the construction of the optimised filter, we 
assume an individual cluster situated at the North pole (dg = 0) 
with a characteristic angular SZ-signal yy(6 = \G\) = ATy(6), where 
we separate the true amplitude A and the spatial profile normalised 
to unity, Ty(9). The underlying cluster profile p(d) is assumed to 
follow a generalised King-profile with an exponent A which is a pa- 
rameter in our analysis. At each observation frequency this profile 
is convolved with the (Gaussian) beam of the respective Planck- 
channel (c.f. Sect. 13. H yielding: 



T,(0) = ^ AQ.'p{e')by(\e-e'\) = Y,Ti^j'; 

f-(\ 



P(0) 



1 + 1^ 



and Tfo., 



47t 

2771' 



(35) 



yP(Q. (36) 



For the second step in eqn. <36t we used the convolution theorem 
on the sphere to be derived in AppendixlAl The background n,,(0) is 
assumed to be a compensated homogeneous and isotropic random 
field with a cross power spectrum Cf,^] w defined by 



(m.yi'^e'm'n/ 



■ Q,y]y2^«'^"!."!" 



where (n„(f )) = 0, 



(37) 



nf„,v denotes the spherical harmonics expansion coefficient of 
ny(0), 6[[> denotes the Kronecker symbol, and (•) corresponds to an 
ensemble average. Assuming ergodicity of the field under consid- 
eration allows taking spatial averages over sufficiently large areas 
CI = <9(4;r) instead of performing the ensemble average. 



6.2 Concepts in filter construction 

The idea of an optimised matche d filter for m ultifrequency observa- 
tions was recently proposed bv lHerranz et aL l2002a) for the case 
of a flat geometry. For each observing frequency, we aim at con- 
structing a centrally symmetric optimised filter function if/y{6) op- 
erating on a sphere. Its functional behaviour induces a family of 
filters if/y(d,Ry) which differ only by a scaling parameter Ry. For 
a particular choice of this parameter, we define the filtered field 
Uy(Ry, f3) to be the convolution of the filter function with the ob- 
served all-sky map at frequency y. 



Uy(Ry,i3) = J dasy(e)4/y(\e- i3\,Ry} 

CO +e 



47r 



'^Cm.y^'eO.viRy) ■ 



(38) 



(39) 



(40) 



For the second step, the convolution theorem to be derived in 
Appendix|^was used. The combined filtered field is defined by 
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(41) 



Taking into account the vanishing expectation value of the noise 
{ny(9)) = 0, the expectation value of the filtered field at the North 
pole /3 = is given by 



0)) 



i/^m.viRy)- 



(42) 



The assumption that the cross power spectrum of the signal is neg- 
ligible compared to the noise power spectrum is justified because 
the thermal and kinetic amplitudes are small compared to unity, 
Ay^„ <c 1. Thus, the variance of the combined filtered field <4U is 
determined by 



■ ,Rn) 



{MR, 



,R^;f3)-{u(Ri,...,R^;f3))f 



= ^^Ce, V, V2 'f'eo, V, (Rvt ) lAfo, V2 (Rvn)- (43) 

.',.■'2 t=o 

The optimised filter functions if/y{d) are chosen to detect the 
clusters at the North pole of the sphere (to which they have been 
translated). They are described by a singly peaked profile which 
is characterised by the scale Rf ^ as given by eqn. <35t . While the 
optimised matched filter is defined to obey the first two of the fol- 
lowing conditions, the optimised scale-adaptive filter is required to 
obey all three conditions: 



.(0). 



(i) The combined filtered field u{R 
estimator of the source amplitude A, i.e. {u{R\ 

(ii) The variance of u(Ri, . . . ,Ri^; (3) has a minimum at the 



p(0) 



0) is an unbiased 
A. 



.,<';0)) 



scales ensuring that the combined filtered field is an 

efficient estimator. 

(iii) The expectation value of the filtered field at the source po- 
sition has an extremum with respect to the the scale Rf \ implying 



.(0) 



i(0) 



{Uy(Ry,0)) = 0. 



(44) 



6.3 Matched filter 

For convenience, we introduce the column vectors ipi = [i/'fo.v]. 
Ff = [fyTfQy], and the inverse C^' of the matrix Cf = [Cf^y^y^]. In 
terms of spherical harmonic expansion coefficients, constraint (i) 
reads 



(45) 



Performing functional variation (with respect to the filter func- 
tion ipt) of cr\iRi, . . . ,Rn) while incorporating the (isoperimetric) 
boundary condition J45> through a Lagrangian multiplier yields the 
spherical matched filter tpf 



ip( = aC[ Ff, where a 



f=0 



(46) 



In any realistic application, the cross power spectrum Cf_,,,„ can be 
computed from observed data provided the cross power spectrum 
of the signal is negligible. The quantities a, F(o, and thus ijjfo can 
be computed in a straightforward manner for a specific frequency 
dependence fy and for a model source profile Ty(9). 



6.4 Scale-adaptive filter on ttie sptiere 

The scale-adaptive filter i/'f satisfying all three conditions is given 
by 



ipf = C-j\aFf + Gf), with Gf = [pf yPyl and 



(47) 



d In Tfo.v 
~dTn7" 



fy[2Tm,y + ({TC0,y-Tl^lQ,)]. (48) 



As motivated in Appendix lA5l the logarithmic derivative of t^q with 
respect to the multipole order ^ is a shorthand notation of the dif- 
ferential quotient which is only valid for f » 1. The quantities el- 
and Py are given by the components 

Q' = (A-')oo, y3v = (A-').o, (49) 
where A is the (1+A')x(l+A') matrix with elements 

^00 = £ P'C^'Fp ^ (f?C7')^ (50) 

f=0 (=0 

A«-|;m,.(C7'F,)^, A„„. (C7')^^, (51) 

In these equations, no summation over the indices is implied. 

6.5 Detection fevel and gain 

As described bv lSanz et alj i200ll) . the concept of constructing an 
optimised filter function for source detection aims at maximising 
the signal-to-noise ratio £)„, 



_ {u{Ru...,Rn-M _ . 
= ~ — — ^ 



(52) 



Computing the dispersion of the unfiltered field on the sphere 
yields the signal-to-noise ratio Z), of a signal on the fluctuating 
background: 



vi,V2 f=0 



(53) 



These considerations allow introducing the gain for comparing the 
signal-to-noise ratios of a peak before and after convolution with a 
filter function: 



(54) 



A, a-„(Ru--.,RN} 

If the noise suppression is successful, the gain g will assume 
values larger than one. If the filters are constructed efficiently, they 
are able to reduce the dispersion (cr,XR\, ■ ■ ■ ,Rn) < o"s) while si- 
multaneously retaining the expectation value of the field 1421 . Due 
to the additional third constraint, the scale-adaptive filter is ex- 
pected to achieve smaller gains compared to the matched filter. 

6.6 Numerical derivation of filter kernels 

For the derivation of suitable filter kernels the source profiles are 
assumed to be generalised King-profiles as described by eqn. <36> 
convolved with the respective Planck-heam superimposed on fluc- 
tuating background given by template (5f,„)y-coefficients. The in- 
version of the matrix Cf (cf. eqns. 01 and Ol) 

can be performed 

using either Gauss- Jordan elimination or LU decomposition, which 
both were found to yield reliable results. In the derivation of the 
scale-adaptive filters, however, it is numerically advantageous to 
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artificially exclude the lower multipoles f < 1 from the calcula- 
tion. Due to the sub-millimetric emission of the Milky Way, the 
lower multipoles are very large. Consequently, the corresponding 
i/'fin-coefficients, ^ < 1, have been set to zero, which is not a seri- 
ous intervention since the filters are designed to amplify structures 
at angular scales well below a degree. For consistency, the multi- 
poles below the quadrupole have been artificially removed in the 
derivation of the matched filters as well. 

In contrast to the P/anc/r-simulation pipeline all numerical cal- 
culations presented here are carried out in terms of fluxes measured 
in Jy and not in antenna temperatures for the following reason: 
Cross-power spectra Cc.v^v2 given in terms of antenna temperatures 
are proportional to (vi • viY^ which results in a suppression of the 
highest frequency channels by a factor of almost 10^ compared to 
the lowest frequency channels. 

Furthermore, by working with fluxes instead of antenna tem- 
peratures, the filters for extracting the SZ-signal show frequency 
dependences which can be understood intuitively. The frequency 
dependence is described by eqns. l29l and l30l The normalisation if 
has been chosen to be 1 arcmin", which corresponds to the weakest 
signals Planck will be able to detect. Because of the smallness of 
the source profiles to be detected, the calculations were carried out 
to multipole orders of t^^^ = 4096, which ensures that the beams 
as well as the source profiles are well described. In the plots in 
Sect. 16.71 the filters depicted are smoothed with a moving average 
window comprising eleven bins for better visualisation. 



6.7 Discussion of filter kernels 

6.7.1 Matched filter 

The spherical harmonics expansion coefficients if/ioy following 
from the matched filter algorithm are depicted in Fig. 1151 for 
four frequencies most relevant to SZ-observation, namely for v = 
100 GHz, V = 143 GHz, v = 217 GHz and v = 353 GHz. As back- 
ground noise components the clean COS data set (left column) and 
the exhaustive GAL data set (right column) are contrasted. The fil- 
ter kernels have been derived for optimised detection of sources 
described by a generalised King-profile with angular core radii 
6c = 3'.0 and 9^ = 5'.0 and asymptotic slope A = 1.0. 

The principle how the matched filter extracts the SZ-signal 
from the maps is explained by Fig. 1151 The SZ-profiles the filter has 
been optimised are small structures at angular scales corresponding 
to multipole moments of ^ ^ 10^. In channels below v = 111 GHz, 
the clusters are observed in absorption and the fluxes are decreased. 
For that reason, the filters have negative amplitudes at these angular 
scales for these specific frequencies. At larger scales, the fluctua- 
tions are suppressed by linear combination of the various channels, 
while the filtering functions show very similar shapes. Optimising 
the filters for detection of core radii of 510 instead of 3!0 result 
in a shift of the negative peak at ^ =^ 10^ to smaller multipole or- 
ders. Instrumental noise which is important at even higher mul- 
tipoles is suppressed by the filter's exponential decline at high { 
above { t 2000. The unwanted CMB fluctuations and all Galactic 
contributions at scales larger than the cluster scale are suppressed 
by weightings with varying sign so that the foregrounds are sub- 
tracted at the stage of forming linear combinations of the (S [„,}y- 
coefficients. 

Furthermore, the contours of the matched filter kernels are 
given in Fig. 1151 as functions of both inverse angular scale £ and 
observing frequency v for differing noise contributions. The fig- 
ures compare filters derived for differing background noise com- 
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Figure 16. Matched filter kernels i//v(d) in real space at SZ-frequencies. i.e. 
for V = 100 GHz (solid line), v = 143 GHz (dashed line), v = 217 GHz 
(dotted line) and v = 353 GHz (dash-dotted line), for a data set including 
the CMB, Galactic foregrounds and instiumental noise (thin lines) and for 
a data set containing all Galactic components in addition to the CMB and 
instrumental noise (thick lines). The filter kernel is optimised for the detec- 
tion of a generalised King-profile with core radius 9^ = 5'.Q and asymptotic 
slope A= 1.0. 



positions. The filters shown serve for the optimised detection of 
generalised King-profiles with core radius 6^ = 15^0 and asymp- 
totic slope A = 1.0. These (rather large) values have been chosen 
for visualisation purposes. For clarity, the contour denoting zero 
values has been omitted due to noisy data. In these figures it is ap- 
parent how the filters combine the frequency information in order 
to achieve a suppression of the unwanted foregrounds: At multipole 
moments of a few hundred, the filters exhibit changes in sign, such 
that the measurements at low frequencies are subtracted from the 
measurements at high frequencies in the linear combination of the 
filtered maps. 

Fig. ll6l illustrates the filter kernels <f/v(d) in real space for the 
same selection of frequencies and background noise components 
as given above. The filter kernels il/y(6) have been synthesised from 
the i/'fo.y-coefficients using the alm2grid-utility of the Planck- 
simulation package. Here, the parameters of the King-profile to be 
detected are (9c,A) = (5'.0, 1.0). The filter kernels are similar in 
shape to Mexican-hat wavelets, but show more than one oscilla- 
tion. Their action on the sky maps is to apply high-pass filtering, 
such that all long-wavelength modes are eliminated. At the cluster 
scale, they implement a linear combination of the sky maps that 
aims at amplifying the SZ-signal: The kernels derived for both the 
V = 100 GHz- and v = 143 GHz-channel exhibit a central depres- 
sion which is used to convert the SZ-signal to positive amplitudes. 
The other two channels resemble simple Gaussian kernels which 
smooth the maps to a common effective angular resolution. At fre- 
quencies of y = 217 GHz and v = 353 GHz the most important 
emission feature is Galactic Dust, which is suppressed by the fil- 
ter's small amplitudes. In this way, the weak SZ-signal is dissected. 

In Fig. 1171 filter kernels derived with both algorithms for point 
sources (i.e. with beam profiles of the respective Planck-chaimels) 
are compared, that have been optimised for the detection of varying 
spectral behaviour of the signal, in this case the thermal SZ-effect, 
the kinetic SZ-elfect and a Planckian thermal emitter with a sur- 
face temperature Tsurfacc of 150 K, such as an asteroid or planet. 
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Figure 15. Upper panel: Spherical harmonics coefficients lAfo.v as derived with the matched filter algorithm, for v = 100 GHz (solid line), v = 143 GHz (dashed 
line), V = 217 GHz (dash-dotted line) and v = 353 GHz (dotted hne). The left and right columns compare filter kernels derived for a data set containing only 
the CMB, the SZ-eft'ects and instrumental noise (left column) and for a data set comprising all Galactic foregrounds in addition (right column). The filter 
kernels are optimised for the detection of generalised King-profiles with core radii 6^ = 3'.0 (thin lines), 0^ = 5'.0 (thick lines) and asymptotic slope A = 1.0. 
Lower panel: Contour plot of the spherical harmonics expansion coefficients i^co.v derived with the matched filter algorithm as a function of both the 
multipole moment order t and Planck's observing frequency v. The filter kernels have been derived for a optimised detection of a generalised King-profile 
with (di:,A) = (15^0, 1.0). The contours are Hnearly spaced in ai'sinh(10^ ■ if/(o,v)- 



The filter kernels depicted correspond to observing frequencies of 
V = 143 GHz and v = 217 GHz. The filters clearly reflect the 
spectral behaviour of the emission laws of the sources one aims at 
detecting: While the filter kernels designed for detecting thermal 
SZ-clusters reflect the peculiar change in sign in the SZ-effect's 
frequency dependence, the other two curves show the behaviour 
to be expected for a Planckian emitter and the kinetic SZ-eff^ect, re- 
spectively. Again, the better angular resolution of the v = 217 GHz- 
channel is apparent by the shifting of the curves to higher multipole 
order £. 



6. 7. 2 Scale-adaptive filter 

The spherical harmonics expansion coefficients i/'fo.v following 
from the scale- adaptive filter algorithm for the frequencies v = 
100 GHz, V = 143 GHz, v = 217 GHz and v = 353 GHz are 
given in the upper panel of Fig. 1181 The left and right columns com- 



pare the filter kernels for differing noise components. Their func- 
tional shape has a number of important features in common with the 
matched filters: They suppress the uncorrected pixel noise, which 
is dominant at high d by their exponential decline at f i 2000. Fur- 
thermore, the filters amplify the SZ-signal, which is negative at fre- 
quencies below y = 217 GHz, by assuming large negative values 
and hence converting the SZ-signal to yield positive amplitudes. 
Additionally, the filters show a distinct secondary peak at ^ ^ 2000 
which causes the kernels to be more compact after transformation 
to real space and enables the size measurement. A more general 
observation is that the scale-adaptive filter kernel shapes are more 
complex and noisier in comparison to the matched filter, especially 
at high {. 

The scale-adaptive filter makes even stronger use of the spec- 
tral information than the matched filter. Especially the contour plots 
in Fig. 1181 show that the scale-adaptive filter exhibits alternating 
signs when varying the observing frequency v while keeping the 
angular scale H fixed. In this way, the noise contributions are iso- 
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Figure 18. Upper panel: Spherical harmonics coefficients i/ro.v as derived with the scale-adaptive filter algorithm, for v = 100 GHz (solid line), v = 143 GHz 
(dashed line), v = 217 GHz (dash-dotted line) and v = 353 GHz (dotted line). The filter kernel is optimised for the detection of generalised King-profiles with 
core radii = 3'0 (thin lines) and 0^ = 5fO (thick lines) and asymptotic slope /i = 1.0. 

Lower panel: Contour plots of the spherical harmonics expansion coefficients i^fo.r derived with the scale-adaptive filter algorithm as a function of both the 
multipole moment order t and Planck's observing frequency v are shown. The filter kernels have been derived for a optimised detection of a generalised 
King-profile with (ff,/l) = (15^0, 1.0). The contours are linearly spaced in arsinh(10^ • ipio^y). Again, the columns compare the influence of the fluctuating 
background, comprising cosmological contributions (CMB and both SZ-eft'ects) and instrumental noise (left column) to the data set that contains all Galactic 
foregrounds in addition. 



lated in angular scale and subsequently suppressed by linear com- 
bination of the maps. Furthermore, one notices a change in sign 
at multipole order C ^ 200 which is common to the frequencies 
V = 100 . . . 353 GHz, at which the CMB signal is strongest. Aim- 
ing at reducing the variance of the filtered maps, the scale-adaptive 
filter is suppressing the (5f„,)y-coefficients by assuming small val- 
ues. 

Fig. 1191 gives the filter kernels 4/y(6) in real space for se- 
lected frequencies and background noise components. The scale- 
adaptive filters work similarly as the matched filters like Mexican- 
hat wavelets and subject the sky maps to high pass filtering. 

In Fig. l20l filter kernels derived with both algorithms for point 
sources (i.e. with beam profiles of the respective Planck-c\iaxmA%) 
are compared, that have been optimised for the detection of varying 
spectral behaviour of the signal, in this case the thermal SZ-effect, 
the kinetic SZ-effect and a Planckian thermal emitter with a surface 
temperature T^adixis of 150 K, such as an asteroid or planet. The 



filter kernels depicted correspond to observing frequencies of y = 
143 GHz and v = 217 GHz. As in the case of the matched filter, 
the frequency dependence of the signal is reflected by the sign of 
the filter kernel at the anticipated angular scale of the profile to be 
detected. 



6.8 Filter renormalisation and synthesis of likelihood maps 

Once the filter kernels are derived, the filtered fields Uy(Ry,J3) can 
be synthesised from the j(f„, ^-coefficients (defined in eqn. <40» and 
the resulting maps can be added in order to yield the co-added, fil- 
tered field ii{Ri ,Rn,P) (see eqn. <4U '). which can be normalised 
by the level of fluctuation (T„ (given by eqn. <42» to yield the like- 
lihood map D(0). It is favourable to divide the filter kernels by the 
variance cr,, and to apply a renormalisation: 
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Figure 17. Comparison of filter kernel i^fo.v-coeflicients derived for dif- 
fering spectral dependences of the signal: thermal SZ-effect (solid line), 
kinetic SZ-efl!'ect (dashed line) and a Planckian emitter with surface tem- 
perature of Tsurface = 150 K (dash-dotted line). All sources are assumed to 
be point-like, i.e. they appear to have the shape of the Planck-heam. The 
curves are given for observing frequencies of y = 143 GHz (thin line) and 
V = 217 GHz (thick line) and have been deiived with the matched filter algo- 
rithm. The noise is a composite of CMB fluctuations. Galactic and ecliptic 
foregrounds and instrumental noise. 
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Figure 20. Comparison of filter kernel i^fo.v-coefficients derived for dif- 
fering spectral dependences of the signal: thermal SZ-effect (solid line), 
kinetic SZ-eft'ect (dashed line) and a Planckian emitter with surface tem- 
perature of Tsurface = 150 K (dash-dotted line). All sources are assumed to 
be point-like, i.e. they appear to have the shape as the Planck-heam. The 
curves are given for observing frequencies of v = 143 GHz (thin line) and 
y = 217 GHz (thick line) and have been derived with the scale-adaptive fil- 
ter algorithm. The noise is a composite of CMB fluctuations. Galactic and 
ecliptic foregrounds and instrumental noise. 
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Figure 19. Scale-adaptive filter kernels i^v(6) in real space, for v = 100 GHz 
(solid line), v = 143 GHz (dashed line), y = 217 GHz (dotted line) and v = 
353 GHz (dash-dotted line), for a data set incorporating the CMB, Galactic 
foregrounds and instrumental noise. The filter kernel is optimised for the 
detection of a generalised King-profile with parameters (9^, /l) = (5^0, 1.0). 



(55) 



In this case, the filter kernels are invariant under changes in pro- 
file normalisation. With these kernels, the filtered flux maps can 
be synthesised from the set of (5f„,)v-coefficients and the resulting 
maps can be co-added to yield the final normalised likelihood map 
D(fi). It is computationally advantageous, however, to interchange 



the last two steps, 

"08) 1 



= > , Uyifi) 

cr„ cr„ ^ 



V f=0 m=-C 



4?r .„ , (Aro.y 
\J rm)v- 



2e+l 



(56) 



■JTiP) (57) 



= ZZ V27TT 



1=0 m=-f 



(58) 



and to derive the -coefficients first, such that the synthesis has to 
be performed only once. Due to the restriction to axially symmet- 
ric kernels, the convolution can be carried out using the alm2map- 
utility rather than totalconvolve. 

Fig. l21l gives a visual impression of the capability of the above 
described filtering schemes: The figure shows a 30° x 30° wide field 
at the ecliptic North pole at a frequency of v = 353 GHz (at the SZ- 
maximum) as observed by Planck, i.e. the image is smoothed to an 
angular resolution of A(? = 5'0 (FWHM) and contains the fluctu- 
ating CMB, all Galactic and ecliptic foregrounds as well as pixel 
noise. Matched and scale-adaptive filter kernels were derived for 
isolating point sources, i.e. for sources that appear to have profiles 
equal to Planck'?, beams of the corresponding channel. For clarity, 
only amplitudes exceeding a threshold value of 1.0 are shown. 

For comparison. Fig. |^ shows the same detail of the input 
thermal SZ-map as well. It is immediately apparent that the ob- 
servation of SZ-clusters without foreground- and noise suppression 
is not possible and that one has to rely on filtering schemes. As a 
comparison with Fig. l21l shows. the filters are clearly able to isolate 
the SZ-clusters and to strongly suppress all spurious noise contri- 
butions. The matched filter, however, shows a slightly better perfor- 
mance and yields more significant peaks due to better background 
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Figure 21. Upper left panel: A 30° X 30° wide field centered on the ecliptic North pole as synthesised from a data set containing the CMB, all Galactic 
and ecliptic foregrounds and instrumental noise for an observing frequency of v = 353 GHz is shown. The amplitudes are proportional to arsinh(r/i(v = 
353 GHz)/|aK) and the field is smoothed with the corresponding Planck-beam of A6 = 5fO (FWHM). Upper right panel: The same field is shown after 
reconstruction from the Df„,-coefficients. Here, filters derived with the matched filter algorithm optimised for detecting point sources were employed. The 
amplitudes are given in detection significances and the shading scales linearly. Lower right panel: Again, the same field is shown after synthesis from the 
Z)ft„-coelficients but in this case, filters derived with the scale-adaptive filter algorithm optimised for detecting point sources were used. The amplitudes are 
stated as detection significances and the shading is linear. In the lower left panel, the corresponding field taken from the oiiginal thermal SZ-map is given for 
comparison. The amplitudes are oc arsinhClO"* ■ y). 



suppression. There are weak residuals present in both maps due 
to incomplete foreground reduction. These residuals however, have 
small amplitudes compared to the SZ-detections. The highest peaks 
exhibit detection significances amounting to 10.6(r in the case of 
the matched filter and 9. Itr in the case of the scale-adaptive filter. 



It should be emphasised that the filters work exceptionally 
well despite the fact that the Milky Way clearly is a non-Gaussian 
feature, whereas Gaussianity of the fluctuating background was an 
important assumption in the derivation of the filter kernels. Fur- 
thermore, the filters sucessfully separate and amplify the weak SZ- 
signal in the presence of seven different noise contributions (CMB, 



Detecting SZ-Clusters with Planck 19 



four Galactic foregrounds, thermal emission from bodies of the 
Solar system and instrumental noise) that exhibit different spec- 
tral behaviours by relying on just nine broad-band measurements. 
Fig. |22| summarises all steps involved in the simulation of Planck- 
observations, filter derivation and signal extraction. 



7 SUMMARY AND CONCLUSION 

A simulation for assessing Planck'?, SZ-capabilities in the pres- 
ence of spurious signals is presented that combines maps of the 
thermal and kinetic SZ-effects with a realisation of the cosmic mi- 
crowave background (CMB), in addition to Galactic foregrounds 
(synchrotron emission, free-free emission, thermal emission from 
dust, CO-line radiation) as well as the sub-millimetric emission 
from celestial bodies of our Solar system. Additionally, observa- 
tional issues such as the finite angular resolution and spatially non- 
uniform instrumental noise of Planck are taken into account. 

• Templates for modelling the free-free emission and the carbon 
monoxide-line emission have been added to the PZancfc-simulation 
pipeline. The free-free template relies on an /f„-survey of the Milky 
Way. The spectral properties of both foregrounds are modelled with 
reasonable parameter choices, i.e. Tp = W* K for the free-free 
plasma temperature and Tco = 20 K for the mean temperature of 
giant molecular clouds. 

• An extensive package for modelling the sub-millimetric emis- 
sion from planet and asteroids has been implemented for Planck, 
that solves the heat balance equation of each celestial body. It takes 
the movement of the planets and asteroids into account, which 
causes, due to Planck's scanning strategy, double detections sep- 
arated by approximate half-year intervals. The total number of as- 
teroids implemented is ^ 1200. 

• The foregrounds have been combined under proper inclusion 
of Planck's frequency response windows in order to yield a set 
of flux maps. The auto- and cross-correlation properties of those 
maps are investigated in detail. Furthermore, their decomposition 
into spherical harmonics <5f„,)y serve as the basis for the filter con- 
struction. It should be emphasised that the spectral properties of a 
foreground component were assumed to be isotropic. 

• In order to separate the SZ-Signal and to suppress the fore- 
ground components, the theory of matched and scale-adaptive fil- 
tering has been extended to spherical data sets. The formulae in the 
context of spherical coordinates and /"'-decomposition are analo- 
gous to those derived for Cartesian coordinate systems and Fourier- 
transforms. 

• The global properties of filter kernel shapes are examined as 
functions of observing channel, composition of noise, parameters 
of the profile to be detected and spectral dependence of the signal. 
Transformation of the filter kernels to real space yields functions 
that resemble the Mexican-hat wavelets, but show more than one 
oscillation. The shape of the filter kernels can be understood intu- 
itively: They subject the maps to high-pass filtering while retaining 
structures similar in angular extent to the predefined profile size. 
The signal is then amplified by linear combination of the maps, 
which again is apparent in the sign of the filter kernels. 

• The functionality of the filtering scheme is verified by apply- 
ing them to simulated observations. It is found that the Galactic 
foregrounds can be suppressed very effectively so that the SZ- 
cluster signals can be retrieved. Comparing the two filters, the 
scale-adaptive filter performs not as good as the matched filter, 
which is in accordance to the findings of Herranz et al. ( 2002a b). 
It should be emphasised that for the derivation of the filter kernels 



nothing but a model profile and all cross-power spectra (in Planck's 
case a total number of 45 independent Cf_viv2"Ssts) are used. 

The scientific exploitation of our simulation and the character- 
isation of Planck's SZ-cluster sample, e.g. the number density as a 
function of detection sigificance as well as filter parameters, spatial 
distribution in dependence on Galactic and ecliptic latitude and the 
distribution in redshift, mass and apparent size will be the subject 
of a forthcoming paper. 
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APPENDIX A: OPTIMISED FILTER FOR SINGLE 
FREQUENCY ALL-SKY OBSERVATIONS 

This appendix derives the optimised filters for single frequency all- 
sky observations and thus serves as a detailed supplement to Sect.|S| 
where optimised filters for multi-frequency observations were de- 
rived. The formalism outlined in this appendix might be applied to 
future all-sky surveys in the X-ray or microwave regime. 



Al Assumptions and definitions 

In order to construct our filters, we consider an all-sky map of the 
detected scalar field s(d) 

s(e) = y(\e - 6>o|) + n(e), (Al) 

where G = (&, ip) denotes a two-dimensional vector on the sphere 
and 00 is the source location. The first term of the right-hand side 
represents the amplitude of the sources to be detected, while the 
second term in corresponds to the generalised noise present in the 
map which is composed of any detected features other than the 
desired signal including for instance instrumental noise. The sta- 
tistical properties of the noise are assumed to be characterised by 
its power spectrum («fm'if/„,/) = Cr5r('6,„j„> ■ In order to sketch the 
construction of the optimised filter, we assume an individual clus- 
ter situated at the North pole {do = 0) with a characteristic angular 
SZ-signal y(9 = \d\) = Arid), separating the amplitude A and the 
profile T(d). 



A2 Convolution theorem on the sphere 

Filtering a scalar field on the sphere with an arbitrary, asymmetric 
kernel requires the specification of the convolution path as well as 
the orientation of the filter kernel at each position on the sphere 
in order to apply any convolution algorithm ( Wandelt & Gors^ 
[2001 ). Because of the simplifying restriction to centrally symmet- 
ric filter kernels, we give the theorem for the convolution of two 
functions, one of which is assumed to be centrally symmetric. The 
filtered field u(f3) is obtained by convolution of the centrally sym- 
metric filter function tf/(d) with the scalar field on the sphere s(d), 

i((/3)= r dn s(eme-f3\). (A2) 
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Expansion of these two scalar fields into spherical harmonics yields 



(A3) 



The last step assumes central symmetry. In this case, only modes 
with m = are contributing. For proceeding, the additio n theorem 
for Legendre polynomials Pc{x) iArfken & WebeJl995h is used in 
substituting y = \6 - fi\: 



Peicosy) 



An 



(A5) 



Combining these equations and applying the completeness relation 
yields the convolution relation for a centrally symmetric filter ker- 
nel. 



= £ £ ^i,x'(fi\ with 



An 



S[„ if/(o ■ (A6) 



2e+ 1 



4n 



-il/m(R)Pr(cose), 



1 f , 2{+\ ie\ 



f(cos ff). 



(A7) 



(A8) 



while allowing for central symmetry of the filter function. For a 
particular choice of R the filtered field u{R, j3) is obtained by con- 
volution (c.f. Appendix. I A2t : 



+[ 



u(R,f3) = 2^ 2^ Ui,„Y';-(f3\ and 

f=0 in=-C 



An 



Stm i//m(R) ■ 



(A9) 



(AlO) 



Taking into account the vanishing expectation value of the noise 
(ny(6)) = 0, the expectation value of the filtered field at the North 
pole /3 = is given by 



{u(RM=AY^Tmikco(R). 



(All) 



A3 Concepts of optimised filtering on the spliere 

The idea of optimised matched filters was proposed by 
iTegmark & de Oliyeira-Cost^ ^^9M), and generalised to scale- 
adaptive filters bv Tsan^^lT j200ll) for a flat topology. The con- 
struction of a centrally symmetric optimised filter function ij/{ff) for 
the amplification and detection of signal profiles differing only in 
size but not in shape implies a family of filters ifi(9/R) introducing 
a scaling parameter R. Decomposing the family of filter functions 
tf/(d/R) into spherical harmonics yields 



Assuming that the power spectrum of the signal is negligible com- 
pared to the noise power spectrum, the variance of the filtered field 
is given by 



a-l(R) = ([u(R,f3) - {u(R,l3))f) = Y. C,,l^l(R). 



(A 12) 



While the optimised matched filter in the case of single fre- 
quency observations is defined to obey the first two of the following 
conditions, the optimised scale-adaptive filter is required to obey 
all three conditions: 
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(i) The filtered field u{R,Q) is an unbiased estimator of the 
source amplitude A at the true source position, i.e. {u(R, 0)) = A. 

(ii) The variance of u(R, (3) has a minimum at the scale R ensur- 
ing that the filtered field is an efficient estimator. 

(iii) The expectation value of the filtered field at the source po- 
sition has an extremum with respect to the the scale R, implying 



— {u{RM = Q- 



A4 Matched filter 



(A 13) 



In order to derive the matched filter, constraint (i) can be rewritten 
yielding 



2 Tm i/'fi 



= 1. 



(A14) 



Performing functional variation (with respect to the filter function 
tj/) of a-l(R) while incorporating the constraint jA14t through a La- 
grangian multiplier yields the spherical matched filter: 



if/fQ = a — , where a 



" 2 



(A15) 



In any realistic application, the power spectrum Cc can be es- 
timated from the observed data provided the power spectrum of the 
signal is negligible. The quantities a, Tiq, and thus the filter kernel 
if/CQ can be straightforwardly computed for any model source profile 

m. 



AS Scale-adaptive filter 

The next step consists of reformulating constraint (iii) in order to 
find a convenient representation for the application of functional 
variation. The expansion coefficient of the family of filter functions 
can be rewritten yielding 



,(R) = ^ j d-e^(^)r?(e) = j ifp^(fi)Y^(R/3), (a16) 

where /? = 6/R. In general, this substitution is not valid, because 
d^6 = sinSd^^di^. In the case of localised source profiles, the angle 
6 is small for non- vanishing values of if/ justifying the approxima- 
tion sin 9x9. The same argument also applies for the boundaries 
of integration. With the aid of eqn. JAl 1> . condition <A13> reads 



dR 



0. 



(A17) 



Using eqn. <A16> . the derivative now acts on the Legendre polyno- 
mial Pf, 



2e+ 1 

4Tr 



TmJ d^/3if/(J3)P'i{cosRj3)j3 sin R/3 = 0. (A18) 



Using the derivative relation of the Legendre polynomials 
iArfken&WebeJl99.'ll. 



P'c(x) = [xPeix) - Pe^iix)], 

1 - 

one obtains 



(A19) 
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. ... . sinR/3 

[cosRfiPeicosRp) - Pm(cosR/3)] = 0. 



(A20) 



In our case, the angle 6 is small for non- vanishing values of i/r justi- 
fying the approximations s'mRp x Rfi and cosR/3 » 1. Substituting 
back, d-/3 = d'd/R^, introducing x = cos 9 = cos Rfj, and inserting 
the inversion of eqn. <A16> . namely 



(A21) 



C'=0 

one arrives at 



4n 



v-i /2f+l /2^' + l 



47T 



In 



I 



dxP(,(x)[Pe(x)- Pm(x)]=0. 



(A22) 



Applying the orthogonality relation for the Legendre polynomials. 



dxPt{x}Pt>{x) = :-—rSte', 

] Zl + i 



(A23) 



and using the small angle approximation in the second term of 
eqn. <A22t with the same argument as given above, yields the fi- 
nal result 



^ (f/m{R)[Tm + {(Tm - rf_i.o)] = 0. 



(A24) 



Replacing the differential quotient with the corresponding deriva- 
tive is a valid approximation for f » 1. Thus, eqn. jA24t can be 
recast in shorthand notation yielding 



2 + 



d In Tfo 



d\n£ 



0. 



(A25) 



This result could have been obtained independently by attaching 
the tangential plane to the North pole and applying Fourier decom- 
position of the filter function ifi and the source profile t. For that 
reason, it is not surprising that the functional form of this condition 
on the sphere agrees with that obtained by Sanz et al. (2001) for a 
flat topology in two dimensions. 

Performing functional variation (with respect to the filter func- 
tion if/) of crliR) while interlacing the constraints <A14> and <A25> 
through a pair of Lagrangian multipliers yields the spherical scale- 
adaptive filter. 



'I'm 



2b + c- (2a + b) 



d In Tfo 



dln^ 



Zl™ - V — ° 



Tfo dTg 

Q~de 



dine 



(A26) 



(A27) 



As before in the case of the matched filter, the power spectrum 
Cc can be derived from observed data provided the power spectrum 
of the signal is negligible. Assuming a model source profile T(d), 
the quantities Tfo, a, b, c, and finally if/fo can be computed in a 
straightforward way. The derivative of Tm with respect to the mul- 
tipole order f is a shorthand notation of the differential quotient in 
eqn. iX^ . 



REFERENCES 

Aghanim N., Hansen S. H., Lagache G., 2004, ArXiv Astro- 
physics e-prints 



22 B. M. Schdfer, C. Pfrommer, R. M. Hell and M. Bartelmann 



Arfken G. B., Weber H. J., 1995, Materials and Manufacturing 

Processes 
Bartelmann M., 2001, A&A, 370, 754 

Bennett C. L., Fixsen D. J., Hinshaw G., Mather J. C, Moseley 
S. H., Wright E. L., Eplee R. E., Gales J., Hewagama T., Isaac- 
man R. B., Shafer R. A., Turpie K., 1994, ApJ, 434, 587 
Bennett C. L., Halpem M., Hinshaw G., et al. 2003, ApJS, 148, 1 
Bennett C. L., Hill R. S., Hinshaw G., Nolta M. R., Odegard N., 
Page L., Spergel D. N., Weiland J. L., Wright E. L., Halpem 
M., Jarosik N., Kogut A., Limon M., Meyer S. S., Tucker G. S., 
WoUack E., 2003, ApJS, 148, 97 
Bersanelli M., Mandolesi N., 2000, Astrophysical Letters Com- 
munications, 37, 171 
Birkinshaw M., 1999, Phys. Rep., 310, 98 
Bouchet F. R., Gispert R., 1999, New Astronomy, 4, 443 
Bouchet R R., Prunet S., Sethi S. K., 1999, MNRAS, 302, 663 
Cay6n L., Sanz J. L., Barreiro R. B., Martinez-Gonzalez E., Vielva 
R, Toffolatti L., Silk J., Diego J. M., Argueso P., 2000, MNRAS, 
315,757 

Colberg J. M., White S. D. M., Yoshida N., MacFarland T. J., 
Jenkins A., Frenk C. S., Pearce F. R., Evrard A. E., Couchman 
H. M. P., Efstathiou G., Peacock J. A., Thomas P. A., The Virgo 
Consortium 2000, MNRAS, 319, 209 

Cremonese G., Marzari F., Burigana C, Maris M., 2002, New As- 
tronomy, 7, 483 

Dame T. M., Hartmann D., Thaddeus P., 1996, Bulletin of the 

American Astronomical Society, 28, 1362 
Dame T. M., Hartmann D., Thaddeus P, 2001, ApJ, 547, 792 
Diego J. M., Vielva P., Martinez-Gonzalez E., Silk J., Sanz J. L., 

2002, MNRAS, 336, 1351 
Dore O., Teyssier R., Bouchet F. R., Vibert D., Prunet S., 2001, 

A&A, 374, 358 
Finkbeiner D. P, 2003, ApJS, 146, 407 

Finkbeiner D. P, Davis M., Schlegel D. J., 1999, ApJ, 524, 867 
Finkbeiner D. P, Davis M., Schlegel D. J., 2000, ApJ, 544, 81 
Geisbiisch J., Kneissl R., Hobson M., 2004, ArXiv Astrophysics 
e-prints 

Giardino G., Banday A. J., Gorski K. M., Bennett K., Jonas J. L., 
Tauber J., 2002, A&A, 387, 82 

Goldin A. B., Kowitt M. S., Cheng E. S., Cottingham D. A., 
Fixsen D. J., Inman C. A., Meyer S. S., Puchalla J. L., Ruhl J. E., 
Silverberg R. R, 1997, ApJL, 488, L161-F 

Griffin M. J., Ade P. A. R., Orton G. S., Robson E. I., Gear W. K., 
Nolt I. G., Radostitz J. V, 1986, Icarus, 65, 244 

Haslam C. G. T., Klein U., Salter C. J., StofFel H., Wilson W. E., 
Cleary M. N., Cooke D. J., Thomasson P, 1981, A&A, 100, 209 

Haslam C. G. T., StofFel H., Salter C. J., Wilson W. E., 1982, A&A 
Suppl. Ser., 47, 1 

Herranz D., Sanz J. L., et al. 2002a, MNRAS, 336, 1057 

Herranz D., Sanz J. L., et al. 2002b, ApJ, 580, 610 

Hivon E., Wandelt B. D., Gorski K. M., 1998, in Evolution of 
Large-Scale Structure: From Recombination to Garching Cos- 
mic Microwave Backgrovmd Anisotropy Power Spectrum Statis- 
tics for High Precision Cosmology 

Hobson M. P., Barreiro R. B., Toffolatti L., Lasenby A. N., Sanz 
J., Jones A. W., Bouchet F R., 1999, MNRAS, 306, 232 

Hobson M. P., Jones A. W., Lasenby A. N., Bouchet F. R., 1998, 
MNRAS, 300, 1 

Jenkins A., Frenk C. S., White S. D. M., Colberg J. M., Cole S., 
Evrard A. E., Couchman H. M. R, Yoshida N., 2001, MNRAS, 
321, 372 

Lagache G., 2003, A&A, 405, 813 



Lagache G., Puget J. L., 2000, A&A, 355, 17 
Lagerros J. S. V, 1996a, A&A, 310, 1011 
Lagerros J. S. V, 1996b, A&A, 315, 625 
Lagerros L S. V, 1997, A&A, 325, 1226 
Lagerros J. S. V, 1998, A&A, 332, 1123 

Lamarre J. M., Puget J. L., Bouchet F., 2003, New Astronomy 

Review, 47, 1017 
Leinert C, Abraham P., Acosta-Pulido J., Lemke D., Siebenmor- 

gen R., 2002, A&A, 393, 1073 
MuUer T. G., 2001, Planetary and Space Science, 49, 787 
MuUer T. G., Lagerros J. S. V., 1998, A&A, 338, 340 
Miiller T. G., Lagerros J. S. V, 1999, AAS/Division for Planetary 

Sciences Meeting, 31,0 
MuUer T. G., Lagerros J. S. V, 2002, A&A, 381, 324 
Maino D., Burigana C, Maltoni M., Wandelt B. D., Gorski K. M., 

Malaspina M., Bersanelli M., Mandolesi N., Banday A. J., Hivon 

E., 1999, A&A Suppl. Ser., 140, 383 
Maino D., Farusi A., Baccigalupi C, Perrotta F, Banday A. J., 

Bedini L., Burigana C, De Zotti G., Gorski K. M., Salerno E., 

2002, MNRAS, 334, 53 
Mather J. C, Fixsen D. J., Shafer R. A., Mosier C, Wilkinson 

D.T, 1999, ApJ, 512, 511 
Moore P., 2000, The data book of astronomy. The data book of 

astronomy, Bristol: Institute of Physics Publishing (lOPP), 2000, 

529 p. ISBN 0750306203 
Natoli R, de Gasperis G., Gheller C, Vittorio N., 2001, A&A, 

372, 346 

Neugebauer G., Miinch G., Kieffer H., Chase S. C, Miner E., 

1971, AJ, 76,719 
Puget J.-L., Abergel A., Bernard J.-P, Boulanger R, Burton W. B., 

Desert F-X., Hartmann D., 1996, A&A, 308, L5+ 
Reach W. T, Morris P., Boulanger F, Okumura K., 2003a, ArXiv 

Astrophysics e-prints 
Reach W. T., Morris P., Boulanger R, Okumura K., 2003b, Icarus, 

164, 384 

RephaeU Y., 1995, ARA&A, 33, 541 

Sanz J. L., Herranz D., Martmez-Gonzalez E., 2001, ApJ, 552, 
484 

Schafer B. M., Pfrommer C, Bartelmann M., 2004, preprint astro- 
ph/040399 

Schlegel D. J., Finkbeiner D. P, Davis M., 1997, Bulletin of the 

American Astronomical Society, 29, 1354 
Schlegel D. J., Finkbeiner D. R, Davis M., 1998, ApJ, 500, 525 
Seljak U., Zaldarriaga M., 1996, ApJ, 469, 437 
Sokasian A., Gawiser E., Smoot G. F, 2001, ApJ, 562, 88 
Song Y., Cooray A., Knox L., Zaldarriaga M., 2003, ApJ, 590, 

664 

Spergel D. N., Verde L., Penis H. V, Komatsu E., Nolta M. R., 
Bennett C. L., Halpem M., Hinshaw G., Jarosik N., Kogut 
A., Limon M., Meyer S. S., Page L., Tucker G. S., Weiland 
J. L., WoUack E., Wright E. L., 2003, ApJ, submitted, astro- 
ph/0302209 

Springel V, Hemquist L., 2002, MNRAS, 333, 649 

Springel V, Yoshida N., White S. D. M., 2001, New Astronomy, 

6, 79 

Stolyarov V, Hobson M. P., Ashdown M. A. J., Lasenby A. N., 

2002, MNRAS, 336, 97 
Sunyaev R. A., Zel'dovich I. B., 1972, Comments Astrophys. 

Space Phys., 4, 173 
Sunyaev R. A., Zel'dovich I. B., 1980, ARA&A, 18, 537 
TegmarkM., de OUveira-Costa A., 1998, ApJL, 500, L83+ 
Tegmark M., Efstathiou G., 1996, MNRAS, 281, 1297 



Detecting SZ- Clusters with Planck 

Toffolatti L., Argueso Gomez R, de Zotti G., Mazzei P., Frances- 

chini A., Danese L., Burigana C, 1998, MNRAS, 297, 117 
Valls-Gabaud D., 1998, PASA, 15, 111 

Vielva P., Barreiro R. B., Hobson M. P., Martmez-Gonzdlez E., 
Lasenby A. N., Sanz J. L., Toffolatti L., 2001, MNRAS, 328, 1 

Vielva P., Martinez-Gonzalez E., Cayon L., Diego J. M., Sanz 
J. L., Toffolatti L., 2001, MNRAS, 326, 181 

Villa P., Sandri M., Bersanelli M., Butler R. C, Mandolesi N., 
Mennella A., Marti-Canales J., Tauber J., 2003, ArXiv Astro- 
physics e-prints 

Wandelt B. D., Gorski K. M., 2001, Phys. Rev. D, 63, 123002 

White M., Hemquist L., Springel V., 2002, ApJ, 579, 16 

White M., Majumdar S., 2003, ArXiv Astrophysics e-prints 

Wright E. L., 1976, ApJ, 210, 250 

Yvon D., Mayet P., 2004, ArXiv Astrophysics e-prints 

This paper has been typeset from a TgX/ ET^ file prepared by the 
author. 



